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Abstract 

This  Semiannual  Technical  Summary^describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS),  for  the  period_l  April  -  30  September ~1990.  Statistics 
are  also  presented  lor  additional  seismic  stations,  which  through  cooperative 
agreements  with  institutions  in  the  host  countries  provide  continuous  data 
to  the  NORSAR  Data  Processing  Center.  fNPBG);>These  stations  comprise 
the  Finnish  Experimental  Seismic  Array  (FINESA),  the  German  Experimen¬ 
tal  Seismic  Array  (GERESS),  and^ffom^heTair^of  1990^ two  3-component 
station  in  Polandf-%siaz~aiid~Stary'Folwark:>’  This  ‘■SemiannuxtljdR.'eport  also 
presents  statistics  from  operation  of  the  Intelligent  Monitoring  System,(IMS). 

The  IMS  has  been  operated  in  an  experimental  mode  using  NORESS  and  1 
ARCESS  data,  and  the  performance  has  been  very  satisfactory.  _ _ / 

'"The  NORSAR  .Detection  .Processing  system  has  been  operated  through¬ 
out  the. period,  jyith  an  average  uptime  of  98.0%(as  compared  to  91.8%  for  the 
previous  reporting  period.  A  total  of  2359  seismic  events  have  been  reported 
in  the  NORSAR.  monthly  seismic  bulletin.  The  performance  of  the  contin¬ 
uous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC 
has  been  Satisfactory.  A  system  for  direct  retrieval  of  NORSAR  waveform 
data  through  an  X.25  connection  has  been  implemented,  and  has  been  tested 
successfully  for  acquiring  such  data  by  AFTAC.  Processing  of  requests  for 
full  NORSAR/NORESS  data  on  magnetic  tapes  has  progressed  according  to 
established  schedules.  ' 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data 
Processing  Center  <^NDPG^of  NORESS, ^ARCESS  and  FINESA  data  have 7?/^ 
been  conducted 'throughout  the  period,  with  an  average  uptime  of  92.9%  for  1) 
NORESS  and  95.6%  for  ARCESS.  The  field  computer  at  the  FINESA  array 
was  damaged  by  a  thunderstorm  in  May,  and  FINESA  was  out  of  operation 
until  the  end  of  August.  The  Intelligent  Monitoring  System  was  installed  at 
NORSAR  in  December  1989  and  has  been  operated  experimentally  since  1 
January  1990.  Results  of  the  IMS  analysis  for  the  reporting  period  are  given. 

There  have  been  no  modifications  made  to  the  NORSAR  data  acquisition 
system.  The  process  of  evaluating  technical  options  for  upgrading  the  array 
is  continuing.  A  test  of  a  full  subarray  acquisition  system  will  be  performed 
during  the  next  reporting  period. 

The  routine  detection  processing  of  NORESS,  ARCESS  and  FINESA  is 
running  satisfactorily  on  each  of  the  arrays’  SUN-3/280  data  acquisition  sys¬ 
tems.  The  routine  processing  of  FINESA  data  at  NORSAR  is  similar  to  what 
is  done  in  Helsinki.  GERESS  data  acquisition  and  detection  processing  has 
commenced  during  the  period. 
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Maintenance  activities  in  the  period  comprise  preventive/corrective  main¬ 
tenance  in  connection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS. 
In  addition,  the  maintenance  center  has  been  involved  with  modification  of 
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equipment  for  FINES  A  and  preparatory  work  in  connection  with  NORESS  I 

HF  instrumentation.  Other  activities  involved  testing  of  the  NORSAR  com-  ~l - L - 


munications  systems. 


j^,We  have  continued  burwork  aimed  at  evaluating  the  stability  of  RMS  Lg 
for  yield  estimation  purposes.^  Using  yield  data  which  have  recently  become 
^available  in  a  Soviet  publication,  we  have  found  that  NORSAR  RMS  Lg  cor- 
/  relates  significantly  better  with  yield  than  does  world-wide  m t.  These  same 
/  yield  data  have  been  used  to  evaluate  the  NORESS  detection  threshold  for 
,  Semipalatinsk  explosions,  and  the  90  per  cent  threshold  in  terms  of  yield  is 
found  to  be  as  low  as  0.1  kt.  This  low  threshold  assumes  that  the  explosion  is 
fully  coupled  and  that  the  background  noise  is  at  a  normal  level. 


An  initial  assessment  of  results  from  real-time  processing  of  GERESS  array 
data  has  been  conducted.  We  have  found  that  a  “generic”  beam  deployment, 
identical  to  that  in  use  at  NORESS  and  ARCESS,  also  works  well  for  the 
GERESS  array.  As  more  experience  is  gained,  we  expect,  however,  that  some 
of  the  special  signal  and  noise  characteristics  at  the  GERESS  site  will  make 
it  desirable  to  modify  some  parameters  in  order  to  optimize  detection  per¬ 
formance.  So  far,  GERESS  detects  a  number  of  earthquake  and  man-made 
explosions  within  500  km  which  are  not  detected  by  other  regional  array  sta¬ 
tions.  It  also  provides  additional  information  which  will  improve  the  location 
N accuracy  for  events  detected  by  multiple  stations  in  the  regional  network. 


>The  generalized  beamforming  technique  developed  by  Ringdal  and  Kvaerna 
is  now  being  applied  to  the  network  of  4  regional  arrays  (NORESS,  ARCESS, 

■  FINESAj  GERESS ).sjnitial  experience  has  confirmed  the  effectiveness  of  this 
/  method  in  automatically  and  rapidly  providing  initial  phase  associations  and 
epicenter  estimates  for  regional  events.  For  a  five-day  test  period,  several 
hundred  regional  events  were  associated,  most  of  which  were  detected  by  one 
array  only  (P  and  Lg  phases).  A  total  of  37  events  had  4  or  more  phase  de¬ 
tections  by  the  regional  network.  Analyst  inspection  revealed  that  36  of  these 
were  properly  associated  and  located,  whereas  one  represented  a  coincidental 
match  of  unrelated  phase  detections.  Further  work  is  now  being  conducted 
to  fine  tune  the  parameters  used  by  this  method,  and  also  to  consider  incor¬ 
porating  a  dynamic  consistency  check  using  the  thresholds  computed  by  the 
threshold  monitoring  technique  to  assess  the  likelihood  of  phase  associations. 


The  concept  of  threshold  monitoring,  introduced  by  Ringdal  and  Kvasrna, 
is  a  method  of  monitoring  the  seismic  amplitude  levels  for  the  purpose  of  using 
this  information  to  assess  the  largest  size  of  events  that  might  go  undetected  by 
a  given  net  work.  >In  an  effort  to  demonstrate  the  capabilities  of  this  threshold 
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.  '  *  ‘  monitoring  concept,  a  preliminary  version  has  been  implemented  into  the  Intel¬ 
ligent  Monitoring  System  (IMS).  Using  recordings  by  the  NORESS,  ARCESS 
and  FINESA  arrays,  the  method  has  been  applied  for  monitoring  continuously 
the  Novaya  Zemlya  test  site  for  a  full  one- week  period.  It  is  demonstrated  that 
the  implementation  of  the  threshold  monitoring  method  in  the  IMS  system  en¬ 
ables  real-time  operation.  The  displays  provided  by  the  threshold  monitor  are 
particularly  valuable  in  pointing  out  time  intervals  of  special  interest,  thus 
aiding  the  analyst  in  his  work.  The  interesting  intervals  can  then  be  exam¬ 
ined  by  different  processing  techniques  to  locate  and  identify  the  events.  In 
the  actual  one-week  monitoring  experiment,  we  have  found  that  the  combined 
three-array  threshold  for  Novaya  Zemlya  is  mj  =  2.2  or  lower  50  %  of  the  time, 
and  mj  =  2.5  or  lower  99  %  of  the  time.  All  instances  where  the  threshold  ex¬ 
ceeds  mj  =  2.5  correspond  to  actual  interfering  seismic  events,  either  regional 
or  teleseismic.  Methods  are  currently  being  studied  on  how  to  further  ana¬ 
lyze  the  data  during  such  interference,  with  the  aim  to  establish  operationally 
reliable  monitoring  at  the  lowest  possible  threshold. 

The  concept  of  a  very  small  “hybrid”  array/3-component  station  has  been 
studied.  We  have  used  the  NORESS  A-ring  subarray  as  a  test  station  for  a 
one-week  on-line  detection  processing  experiment.  We  have  found  that  this 
small  array,  which  has  a  diameter  of  only  300  m  (i.e.,  1/10  of  that  of  NORESS) 
has  a  remarkably  good  performance  in  distinguishing  P  and  S  phases  (based 
on  phase  velocity  estimates)  and  in  obtaining  useful  azimuth  estimates  for  all 
phase  types.  Thus  it  would  appear  that  supplementing  3-component  stations 
with  small  triangular  arrays  of  the  A-ring  type  would  to  a  large  extent  alleviate 
the  problems  now  encountered  in  3-component  analysis  of  secondary  phases. 

Work  is  proceeding  toward  establishing  two  modern  3-component  stations 
in  Poland  (Ksiaz  and  Stary  Folwark)  with  continuous  data  transmission  by 
satellite  to  NDPC  for  integration  into  the  IMS  processing.  An  initial  report  on 
the  system  design,  including  site  descriptions,  communication  arrangements 
and  the  NDPC  data  acquisition  system  is  presented.  An  example  of  data 
recorded  at  Ksiaz  indicates  that  data  from  this  station  are  of  high  quality  and 
will  give  a  valuable  contribution  to  the  regional  station  network. 

Further  work  on  statistically  optimal  array  detection,  using  a  generaliza¬ 
tion  of  Capon’s  maximum  likelihood  technique,  is  reported.  Examples  are 
presented  both  for  array  detection  and  onset  time  estimation  of  regional  and 
teleseismic  phases.  It  is  shown  in  particular  that  the  very  strong,  low  frequent 
microseismic  noise  observed  at  NORESS  and  ARCESS  can  be  effectively  sup¬ 
pressed,  due  to  the  high  noise  coherency.  Thus  it  is  possible  to  regenerate 
signal  waveforms  for  which  the  signal  is  represented  in  an  undistorted  way, 
whereas  the  microseismic  noise  is  suppressed  by  as  much  as  12-18  dB. 
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1  Summary 


This  Semiannual  Technical  Summary  describes  the  operation,  maintenance 
and  research  activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Nor¬ 
wegian  Regional  Seismic  Array  (NORESS)  and  the  Arctic  Regional  Seismic 
Array  (ARCESS)  for  the  period  1  April  -  30  September  1990.  Statistics 
are  also  presented  for  additional  seismic  stations,  which  through  cooperative 
agreements  with  institutions  in  the  host  countries  provide  continuous  data 
to  the  NORSAR  Data  Processing  Center  (NPDC).  These  stations  comprise 
the  Finnish  Experimental  Seismic  Array  (FINESA),  the  German  Experimen¬ 
tal  Seismic  Array  (GERESS),  and  from  the  fall  of  1990,  two  3-component 
station  in  Poland:  Ksiaz  and  Stary  Folwark.  This  Semiannual  Report  also 
presents  statistics  from  operation  of  the  Intelligent  Monitoring  System  (IMS). 
The  IMS  has  been  operated  in  an  experimental  mode  using  NORESS  and 
ARCESS  data,  and  the  performance  has  been  very  satisfactory. 

The  NORSAR  Detection  Processing  system  has  been  operated  through¬ 
out  the  period  with  an  average  uptime  of  98.0%  els  compared  to  91.8%  for  the 
previous  reporting  period.  A  total  of  2359  seismic  events  have  been  reported 
in  the  NORSAR  monthly  seismic  bulletin.  The  performance  of  the  contin¬ 
uous  alarm  system  and  the  automatic  bulletin  transfer  by  telex  to  AFTAC 
has  been  satisfactory.  A  system  for  direct  retrieval  of  NORSAR  waveform 
data  through  an  X.25  connection  has  been  implemented,  and  has  been  tested 
successfully  for  acquiring  such  data  by  AFTAC.  Processing  of  requests  for 
full  NORSAR/NORESS  data  on  magnetic  tapes  Iiels  progressed  according  to 
established  schedules. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data 
Processing  Center  (NDPC)  of  NORESS,  ARCESS  and  FINESA  data  have 
been  conducted  throughout  the  period,  with  an  average  uptime  of  92.9%  for 
NORESS  and  95.6%  for  ARCESS.  The  field  computer  at  the  FINESA  array 
was  damaged  during  a  thunderstorm  in  May,  and  FINESA  was  out  of  operation 
until  the  end  of  August.  The  Intelligent  Monitoring  System  wels  installed  at 
NORSAR  in  December  1989  and  has  been  operated  experimentally  since  1 
January  1990.  Results  of  the  IMS  analysis  for  the  reporting  period  are  given. 

There  have  been  no  modifications  made  to  the  NORSAR  data  acquisition 
system.  The  process  of  evaluating  technical  options  for  upgrading  the  a: ray 
is  continuing.  A  test  of  a  full  subarray  acquisition  system  will  be  performed 
during  the  next  reporting  period. 

The  routine  detection  processing  of  NORESS,  ARCESS  and  FINESA  is 
running  satisfactorily  on  each  of  the  arrays’  SUN-3/280  data  acquisition  sys¬ 
tems.  The  routine  processing  of  FINESA  data  at  NORSAR  is  similar  to  what 
is  done  in  Helsinki.  GERESS  data  acquisition  and  detection  processing  has 
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commenced  during  the  period. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  main¬ 
tenance  in  connection  with  ell  the  NORSAR  subarrays,  NORESS  and  ARCESS. 
In  addition,  the  maintenance  center  has  been  involved  with  modification  of 
equipment  for  FINESA  and  preparatory  work  in  connection  with  NORESS 
HF  instrumentation.  Other  activities  involved  testing  of  the  NORSAR  com¬ 
munications  systems. 

The  research  activity  is  summarized  in  Section  7.  Section  7.1  presents  a 
continuation  of  the  evaluation  of  the  stability  of  RMS  Lg  for  yield  estimation 
purposes.  Initial  results  from  real-time  processing  of  GERESS  array  data  are 
given  in  Section  7.2.  In  Section  7.3  a  study  of  generalized  beamforming  using  a 
network  of  regional  arrays  is  presented.  The  continuous  threshold  monitoring 
technique  as  applied  to  the  Novaya  Zemlya  test  site  is  discussed  in  Section 
7.4.  Section  7.5  discusses  real-time  processing  using  data  from  a  hybrid  3- 
component/small  array  station.  In  Section  7.6  development  work,  including 
system  design,  site  descriptions,  communication  arrangements  and  the  NDPC 
data  acquisition  system,  in  connection  with  two  new  three-component  stations 
in  Poland  is  presented.  Optimal  group  filtering  and  noise  attenuation  for  the 
NORESS  and  ARCESS  arrays  is  discussed  in  Section  7.7. 
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NORSAR  Operation 


2.1  Detection  Processor  (DP)  operation 

There  have  been  55  breaks  in  the  otherwise  continuous  operation  of  the  NOR¬ 
SAR  online  system  within  the  current  6-month  reporting  interval.  The  uptime 
percentage  for  the  period  is  98.0  as  compared  to  91.8  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  down¬ 
time  for  the  days  between  1  April  and  30  September  1990.  The  monthly 
recording  times  and  percentages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 

a)  Hardware  failure  12 

b)  Stops  related  to  program  work  or  error  4 

c)  Hardware  maintenance  stops  2 

d)  Power  jumps  and  breaks  3 

e)  TOD  error  correction  10 

f)  Communication  lines  24 


The  total  downtime  for  the  period  was  87  hours  and  11  minutes.  The 
mean-time-between-failures  (MTBF)  was  3.2  days,  as  compared  to  1.6  for  the 
previous  period. 


J.  Torstveit 
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Fig.  2.1.1  Detection  Processor  downtime  for  April  (top),  May  (middle)  and 
June  (bottom)  1990. 
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Fig.  2.1.1  Detection  Processor  downtime  for  July  (top),  August  (middle)  and 
September  (bottom)  1990. 
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LIST 

OF  BREAKS 

IN  DP 

DAY 

START 

STOP 

92 

5 

40 

5 

93 

4 

57 

5 

95 

7 

59 

8 

107 

6 

0 

6 

113 

6 

14 

6 

116 

12 

35 

12 

122 

12 

0 

12 

154 

6 

14 

8 

156 

6 

54 

6 

158 

6 

10 

6 

158 

11 

24 

11 

159 

8 

45 

8 

159 

9 

2 

9 

159 

9 

37 

9 

159 

9 

58 

10 

159 

12 

12 

12 

160 

18 

21 

21 

166 

9 

16 

9 

171 

13 

0 

24 

172 

0 

0 

7 

176 

12 

8 

12 

178 

6 

31 

6 

184 

5 

42 

5 

184 

10 

53 

10 

185 

10 

28 

10 

185 

10 

34 

10 

187 

11 

29 

11 

192 

6 

28 

11 

196 

16 

35 

19 

199 

11 

48 

11 

200 

12 

25 

12 

201 

10 

18 

10 

206 

8 

30 

9 

225 

6 

48 

6 

225 

11 

55 

11 

227 

11 

2 

12 

228 

10 

7 

11 

228 

11 

32 

11 

232 

6 

23 

6 

239 

10 

45 

10 

240 

1 

30 

4 

240 

6 

27 

6 

240 

6 

34 

7 

249 

6 

0 

6 

251 

12 

43 

24 

252 

0 

0 

10 

259 

5 

49 

12 

260 

2 

49 

6 

261 

2 

45 

3 

261 

5 

38 

5 

263 

0 

18 

5 

263 

13 

43 

18 

263 

19 

15 

19 

264 

12 

1 

12 

264 

21 

59 

23 

269 

19 

25 

20 

271 

13 

18 

13 

PROCESSING  THE  LAST 
COMMENTS . 

42  TOD  RETARED  14MS 
0  TOD  RETARED  7MS 

6  TOD  ADJUSTMENT 

2  TOD  RETARED  2 7 MS 
16  TOD  RETARED  13MS 

37  TOD  RETARED  40MS 
2  TOD  RETARED  12  S 

19  CPU  FAILURE 
57  TOD  RETARED  7MS 
14  LINE  FAILURE 

28  LINE  FAILURE 
48  LINE  FAILURE 

7  LINE  FAILURE 
40  LINE  FAILURE 

0  LINE  FAILURE 
14  LINE  FAILURE 
36  AIRCOND.  FAILURE 

19  LINE  FAILURE 
0  LINE  FAILURE 

44  LINE  FAILURE 
22  CE  MAINTENANC 

34  TOD  RETARED  3 IMS 

43  LINE  FAILURE 
54  LINE  FAILURE 

29  LINE  FAILURE 

35  LINE  FAILURE 

30  LINE  FAILURE 

27  MODCOMP  FAILURE 
42  AIRCOND.  FAILURE 
50  LINE  FAILURE 
27  LINE  FAILURE 

20  LINE  FAILURE 

19  POWER  FAILURE 
50  LINE  FAILURE 
57  LINE  FAILURE 
47  DP-ONLINE  TEST 

20  DP-ONLINE  TEST 
34  LINE  FAILURE 
25  LINE  FAILURE 
47  LINE  FAILURE 
40  POWER  BREAKE 
30  LINE  FAILURE 

38  POWER  BREAKE 

2  TOD  RETARED  8MS 
0  SOFTWARE  PROBLEM 
6  SOFTWARE  PROBLEM 
25  HARDWARE  FAILURE 

36  HARDWARE  FAILURE 

45  HARDWARE  FAILURE 
40  LINE  FAILURE 

20  SOFTWARE  PROBLEM 
0  MAINTENANCE  POWER 
45  HARDWARE  FAILURE 
15  HARDWARE  FAILURE 
22  HARDWARE  FAILURE 
33  HARDWARE  FAILURE 
25  HARDWARE  FAILURE 


HALF-YEAR 


Table  2.1.1  Daily  DP  downtime  in  the  period  1  April  -  30  September  1990. 
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Month 

DP  Uptime 
Hours 

DP  Uptime 
% 

No.  of 

DP  Breaks 

No.  of  Days 
with  Breaks 

DP  MTBF* 
(days) 

APR  90 

719.72 

99.96 

6 

6 

4.3 

MAY  90 

743.97 

100.00 

1 

1 

JUN  90 

695.17 

96.55 

14 

10 

1.9 

JUL  90 

743.92 

98.88 

11 

9 

2.6 

AUG  90 

736.58 

99.00 

10 

6 

2.8 

SEP  90 

674.50 

93.68 

13 

10 

2.0 

98.03 

55 

42 

3.2 

*Mean-time-between-failures  =  total  uptime/no.  of  up  intervals. 


Table  2.1.2  Online  system  performance,  1  April  -  30  September  1990. 
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2.2 


Array  communications 


General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  through¬ 
out  the  reporting  period. 

April  and  May  were  good  months  regarding  the  communications  systems. 
02C  failed  due  to  loss  of  power,  and  06C  had  problems  due  to  an  irregularity 
related  to  the  CTV  equipment. 

June  was  characterized  by  unstable  line  conditions  affecting  all  the  sub¬ 
arrays,  but  also  subarray  electronics  failed  (02B)  and  a  communication  cable 
between  Berger  central  and  the  CTV  (04C)  failed. 

In  July  02B  lost  power  (week  28),  03C  lost  synchronization  frequently 
(week  28),  and  06C  was  down  approximately  8  hours  (week  27)  in  connection 
with  loss  of  synchronization  after  TOD  adjustment. 

In  August  4  subarrays  were  affected:  02B  (weeks  33,34)  caused  by  lightning 
and  loss  of  power;  02C  (week  33)  due  to  an  unknown  cause;  04C  also  in 
connection  with  lightning  and  loss  of  power;  06C  was  interrupted  11  August, 
also  due  to  an  unknown  reason  and  resumed  operation  13  August  after  a 
Modcomp  restart. 

September  was  also  a  good  period  with  regard  to  the  communications 
systems.  Only  02B  was  affected,  experiencing  spikes  on  SP  channel  04.  As 
previously  stated  in  our  reports,  the  TOD  (Time  of  Day  Generator)  is  less 
reliable  than  before,  and  will  probably  soon  be  replaced  with  by  a  new  system. 

Detailed  summary 

April  (weeks  lj-17),  2-29.4.90 

02C  was  down  from  28  April  due  to  loss  of  power. 

May  (weeks  18-22),  30-4-3.6.90 

As  stated  above  02C  was  affected  by  loss  of  power  (week  22)  and  06C  was 
affected  by  an  irregularity  related  to  the  CTV  equipment  (weeks  20,  21  and 
22). 

June  (weeks  23-26),  4-6-1.7.90 

In  June  all  the  subarray  communications  systems  were  affected  by  the  NTA 
transmission  system  (20-21  June).  In  addition,  individual  subarrays  were  also 
affected:  02B  by  power  outage  and  STEM  failure  (weeks  23,  24  and  25)  and 
02C  (weeks  23,  24)  by  power  failure.  Also  04C  continued  to  fail  in  connection 
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with  a  bad  cable  section  between  Berger  central  and  the  CTV. 

My  (weeks  27-30),  2-29.7.90 

02B,  which  is  in  an  exposed  position  with  respect  to  wind,  lightning  and 
snow  conditions,  again  failed  (9-16  July)  due  to  loss  of  power.  Afterwards 
NMC  Hamar  had  to  reset  the  SLEM.  03C  was  affected  4,  10-11  July  by 
excessive  resynchronization  attempts. 

The  04C  communications  problems,  where  NTA,  NMC  and  NDPC  were 
involved,  were  brought  to  an  end  10  July.  The  next  two  weeks  of  the  period 
the  communications  system  operated  well.  06C  was  affected  4  July  for  ap¬ 
proximately  8  hours  as  a  cable  between  Lpten  and  the  subarray  was  cut  by  an 
excavator. 

August  (weeks  31-35),  30.7-2.9.90 

02B  was  affected  15  and  17  August  by  lightning  and  loss  of  power.  In  both 
cases  visits  to  the  subarray  were  necessary  to  restore  the  SLEM  operation. 

18  and  20  August  02C  went  down,  reason  not  found.  After  Modcomp 
restart,  the  system  resumed  operation. 

Also  04 C  was  affected  by  lightning  and  power  loss  9  and  16  August. 

06C  was  interrupted  11  August  due  to  an  unknown  cause.  After  a  Mod¬ 
comp  restart  13  August,  the  system  resumed  operation. 

September  (weeks  36-39),  3-30.9.90 

In  September  all  communications  systems  functioned  satisfactorily.  02B 
was  affected  between  22  and  24  September  probably  under  the  influence  of 
spikes  on  SP  channel  4. 


O.A.  Hansen 
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Sub- 

Apr  90  (4) 

May  90  (5) 

Jun  90  (4) 

Jul  90  (4) 

Aug  90  (5) 

Sep  90  (4) 

Average 

arrays 

2-29.4 

30.4-3.6 

4.6-1.7 

2-29.7 

30.7-2.9 

3-30.9 

1/2  year 

01A 

0.0020 

0.0030 

*1*0.002 

0.002 

0.012 

0.003 

0.004 

01B 

0.0060 

0.0050 

3>*0.004 

0.009 

0.003 

0.012 

0.006 

02B 

0.0020 

0.0040 

4>*n/a 

9>*0.003 

13>*0.002 

17)*0.014 

18>0.005 

02C 

0.0009 

0.0020 

5>*n/a 

0.006 

l3)*0.007 

0.003 

19)0.004 

03C 

0.0010 

0.0020 

®)*0.600 

10)*0.011 

0.004 

0.022 

0.107 

04C 

0.0010 

0.0010 

7>*0.890 

n)*0.008 

1S))*0.006 

0.001 

0.150 

06C 

0.0010 

^O.OSIO 

3)*0.890 

1J>*0.0007 

16>*0.001 

0.006 

*0.103 

AVER 

0.0020 

0.0070 

0.470 

0.005 

0.005 

0.009 

0.054 

LESS 

02B 

02C 

*See  Section  2.2  regarding  figures  preceded  by  an  asterisk. 

Figures  representing  error  rate  (in  per  cent)  preceded  by  a  number  1),  2),  etc.,  are  related  to 
legend  below. 


1), 6), 11), 13) 

Average  2  weeks  (18,19  /  24,26  /  29,30  /  31,32) 

4), 5), 7), 8) 

Average  1  week  (26,24) 

2), 3), 9), 10), 12) 

Average  3  weeks  (23,24,26  /  27,29,30  /  28-30) 

H),15),16) 

Average  3  weeks  (31,32,35  /  31,34,35) 

17) 

Average  2  weeks  (38,39) 

18),  19) 

Average  5  months 

Table  2.2.1  Communications  performance.  The  numbers  represent  error  rates 
in  per  cent  based  on  total  transmitted  frames/week  (1  April  -  30  September 
1990). 
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2.3  Event  detection  operation 

In  Table  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor 
operation  are  given.  The  table  lists  the  total  number  of  detections  (DPX) 
triggered  by  the  on-line  detector,  the  total  number  of  detections  processed  by 
the  automatic  event  processor  (EPX)  and  the  total  number  of  events  accepted 
after  analyst  review  (teleseismic  phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

Accepted  events 
P-phases  Core  phases 

Sum 

Daily 

Apr  90 

10350 

1282 

300 

70 

370 

12.3 

May  90 

7100 

1209 

348 

87 

435 

14.0 

Jun  90 

7700 

1326 

284 

70 

354 

11.8 

Jul  90 

8100 

1401 

482 

69 

551 

17.8 

Aug  90 

8025 

1203 

287 

79 

366 

11.8 

Sep  90 

8475 

1007 

233 

50 

283 

9.4 

1934 

425 

2359 

12.9 

Table  2.3.1.  Detection  and  Event  Processor  statistics,  1  April  -  30  September 
1990. 


B.  Paulsen 
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3  Operation  of  Regional  Arrays 

3.1  Recording  of  NORGSS  data  at  NDPC,  Kjeiier 

Table  3.1.1  lists  the  main  outage  times  and  reasons,  and  as  can  be  seen  the 
main  reasons  for  the  outages  are  hardware  failure  at  the  HUB  and  line  failure. 

The  average  recording  time  was  92.9%  as  compared  to  96.9%  for  the  pre¬ 
vious  period. 


Date 

Time 

Cause 

3  Apr 

0216-0401 

Power  failure  HUB 

3  Apr 

0640-0808 

Power  failure  HUB 

3  Apr 

0829-1210 

Power  failure  HUB 

3  Apr 

2041-2155 

Power  failure  HUB 

3  Apr 

2159- 

Power  failure  HUB 

4  Apr 

-1107 

Power  failure  HUB 

6  Apr 

2338- 

Software  failure 

7  Apr 

-1415 

Software  failure 

7  Apr 

1742-1752 

Transmission  line  failure 

10  Apr 

0511-0609 

Transmission  line  failure 

10  Apr 

1200-2008 

Hardware  failure  HUB 

21  Apr 

0014-0029 

Transmission  line  failure 

26  Apr 

1937-1959 

Transmission  line  maintenance 

2  May 

1931- 

Hardware  failure  HUB 

8  May 

-1014 

Hardware  failure  HUB 

14  May 

0734- 

Hardware  failure  HUB 

15  May 

-0921 

Hardware  failure  HUB 

29  May 

0039-0537 

Transmission  line  failure 

29  May 

0555-0604 

Transmission  line  failure 

31  May 

1423-1918 

Power  failure  HUB 

5  Jun 

0523-1510 

Software  work 

6  Jun 

1256-1336 

Transmission  line  failure 

11  Jun 

0600-0643 

Transmission  line  failure 

11  Jun 

1106-1339 

Transmission  line  failure 

12  Jun 

2207-2304 

Power  failure  HUB 

14  Jun 

1050-1107 

Hardware  failure  NDPC 

19  Jun 

2207-2319 

Power  break  HUB 

28  Jun 

0859-1506 

Power  break  HUB 

28  Jun 

1728-1753 

Power  break  HUB 
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4  Jul 

0727- 

Hardware  failure  NDPC 

5  Jul 

“  -1718 

Hardware  failure  NDPC 

7  Jul 

1102-1108 

Transmission  line  failure 

25  Jul 

0830-0943 

Power  break  NDPC 

27  Jul 

0004-0057 

Hardware  failure  NDPC 

27  Jul 

0949-1158 

Hardware  failure  NDPC 

1  Aug 

1150-1223 

Hardware  maintenance  HUB 

8  Aug 

1542-1608 

Hardware  failure  NDPC 

17  Aug 

1211-1~35 

Transmission  line  test 

17  Aug 

1259-3334 

Transmission  line  test 

21  Aug 

0829-0902 

Transmission  line  test 

28  Aug 

0156-0302 

Power  break  NDPC 

28  Aug 

0823-1031 

Power  break  NDFC 

3  Sep 

0844-0915 

Power  break  NDPC 

3  Sep 

0952-1010 

Power  break  NDPC 

3  Sep 

1023-1027 

Power  break  NDPC 

13  Sep 

1101-1204 

Transmission  line  test 

17  Sep 

1213-1234 

Transmission  line  test 

26  Sep 

1239-1340 

Hardware  work  HUB 

30  Sep 

0100-0200 

Software  failure 

Table  3.1.1.  Interruptions  in  recording  of  NOEESS  data  at  NDPC,  1  April 
-  30  September  1990. 

Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 


April 

93.6% 

May 

76.5% 

June 

96.3% 

'  July 

93.6% 

August 

98.5% 

September 

98.9% 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 

J.  'forstveit 
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Fig.  3.1.1.  (coat.)  NORESS  data  recording  uptime  for  July  (top),  August 
(middle)  and  September  (bottom)  1990. 
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3.2  Recording  of  ARCESS  data  at  NDPC,  Kjelier 


The  main  reasons  causing  most  of  the  ARCESS  outage  in  the  period  are: 
Hardware  failure  at  NDPC  or  HUB,  power  failure  at  HUB  and  work  by  the 
power  company  to  reinforce  the  line  from  which  power  is  supplied  for  ARCESS. 
Outage  intervals  are  listed  in  Table  3.2.1. 

The  average  recording  time  was  95.6%  as  compared  to  95.3%  for  the  pre¬ 
vious  period. 


Date 

Time 

Cause 

24  Apr 

2058- 

Power  break  HUB 

25  Apr 

-0355 

Power  break  HUB 

25  Apr 

1322-1336 

Power  break  HUB 

25  Apr 

2100- 

Power  break  HUB 

26  Apr 

-0408 

Power  break  HUB 

26  Apr 

2110-0648 

Power  break  HUB 

9  May 

1124-1348 

Software  maintenance 

10  May 

1426- 

Hardware  failure  NDPC 

11  May 

-0451 

Hardware  failure  NDPC 

11  May 

0641-0819 

Hardware  failure  NDPC 

13  May 

2102- 

Work  on  powerline  HUB 

14  May 

-0605 

Work  on  powerline  HUB 

14  May 

2101- 

Work  on  powerline  HUB 

15  May 

-0527 

Work  on  powerline  HUB 

15  May 

2100- 

Work  on  powerline  HUB 

16  May 

-0535 

Work  on  powerline  HUB 

20  May 

2102- 

Work  on  powerline  HUB 

21  May 

-0324 

Work  on  powerline  HUB 

21  May 

2102- 

Work  on  powerline  HUB 

22  May 

-0608 

Work  on  powerline  HUB 

22  May 

2058- 

Work  on  powerline  HUB 

23  May 

-0451 

Work  on  powerline  HUB 

25  May 

2219-2319 

Hardware  failure  NDPC 

26  May 

1315-1348 

Hardware  failure  NDPC 

6  Jun 

1747-1916 

Transmission  line  failure 

7  Jun 

1540-1724 

Transmission  line  failure 

8  Jun 

1019-1052 

Hardware  failure  NDPC 

10  Jun 

1149-1210 

Hardware  failure  NDPC 

1  Jul 

0828-0912 

Software  failure 

3  Jul 

1437- 

Transmission  line  failure 
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4 

Jul 

-1011 

Transmission  line 

i  failure 

5 

Jul 

1058-1115 

Software 

failure 

5 

Jul 

1502-1553 

Software 

failure 

7 

Jul 

0850-0913 

Software 

failure 

7 

Jul 

1434-1447 

Hardware 

failure 

NDPC 

8 

Jul 

0902-0925 

Hardware 

failure 

NDPC 

9 

Jul 

2112-2125 

Software 

failure 

14 

Jul 

0504-0620 

Software 

failure 

14 

Jul 

1730-1738 

Software 

failure 

21 

Jul 

2134- 

Software 

failure 

22 

Jul 

-0809 

Software 

failure 

23 

Jul 

1156- 

Hardware 

failure 

HUB 

24 

Jul 

-1854 

Hardware 

failure 

HUB 

25 

Jul 

0757-0839 

Power  break  NDPC 

3 

Aug 

0705-0814 

Software 

failure 

4 

Aug 

0942-0952 

Software 

failure 

22 

Aug 

0437-0522 

Software 

failure 

23 

Aug 

0531-0556 

Software 

failure 

23 

Aug 

1832-1849 

Software 

failure 

28 

Aug 

0114-0305 

Power  break  NDPC 

28 

Aug 

0548-0715 

Power  bread?  NDPC 

28 

Aug 

0823-1031 

Power  break  NDPC 

3 

Sep 

0844-0915 

Power  break  NDPC 

3 

Sep 

0952-1010 

Power  break  NDPC 

3 

Sep 

1023-1027 

Power  bread?  NDPC 

13 

Sep 

2238- 

Hardware 

failure 

NDPC 

14 

Sep 

-0529 

Hardware 

failure 

NDPC 

20 

Sep 

1329-1339 

Work  on  power  NDPC 

20 

Sep 

1632-1749 

Power  break  HUB 

25 

Sep 

1529- 

Hardware 

failure 

NDPC 

26 

Sep 

-0539 

Hardware  failure 

NDPC 

Table  3.2.1.  The  main  interruptions  in  recording  of  ARCESS  data  at  NDPC, 
1  April  -  30  September  1990. 

Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 
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\ 

I 


April 

96.7% 

May 

90.7% 

June 

99.4% 

July 

91.1% 

August 

98.9% 

September 

96.7% 

Fig.  3.2.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 


J.  Torstveit 
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Fig.  3.2.1.  ARCESS  data  recording  uptime  for  April  (top),  May  (middle) 
and  June  (bottom)  1990. 
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3.3 


Recording  of  FINESA  data  at  NDPC,  Kjeller 


The  average  recording  time  was  39.0%  as  compared  to  98.5%  for  the  previous 
period.  The  main  reason  for  the  low  uptime  of  the  FINESA  array  was  a 
thunderstorm  on  May  22  that  put  the  field  computer  out  of  operation  until 
August  26  except  for  a  few  hours  on  June  19. 


Date 

Time 

Cause 

2  May 

0625- 

Hardw are/software  problems  NDPC 

8  May 

-1925 

Hardware/software  problems  NDPC 

15  May 

0044-0119 

Transmission  line  failure 

15  May 

1653- 

Transmission  line  failure 

16  May 

-1103 

Transmission  line  failure 

22  May 

0926- 

Hardware  failure  HUB 

19  Jun 

-0752 

Hardware  failure  HUB 

19  Jun 

1231- 

Hardware  failure  HUB 

26  Aug 

-0830 

Hardware  failure  HUB 

26  Aug 

2328- 

Hardware  failure  HUB 

29  Aug 

-0955 

Hardware  failure  HUB 

29  Aug 

1002-1009 

Hardware  failure  HUB 

31  Aug 

1132-1956 

Hardware  failure  HUB 

6  Sep 

0940- 

Transmission  line  failure 

12  Sep 

-0805 

Transmission  line  failure 

13  Sep 

2314-2332 

Transmission  line  failure 

14  Sep 

0030-0054 

Transmission  line  failure 

18  Sep 

0515-0520 

Transmission  line  failure 

20  Sep 

1328-1400 

Maintenance  NDPC 

20  Sep 

1718-1737 

Maintenance  NDPC 

26  Sep 

0958-1019 

Transmission  line  failure 

26  Sep 

1037-1154 

Transmission  line  failure 

27  Sep 

0743-1302 

Transmission  line  failure 

Table  3.3,1  The  main  interruptions  in  recording  of  FTNESA  data  at  NDPC, 
1  April  -  30  September  1990. 

Monthly  uptimes  for  the  FINESA  on-line  data  recording  task,  taking  into 
account  all  factors  (field  installations,  transmissions  line,  data  center  opera¬ 
tion)  affecting  this  task  were  as  follows: 
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July 

August 

September 


0.0% 

9.2% 

79.0% 


Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently, 
the  availability  of  FINESA  data  in  our  tape  archive,  on  a  day-by-day  basis, 
for  the  reporting  period. 


J.  Torstveit 
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Fig.  3.3.1.  FINESA  data  recording  uptime  for  April  (top),  May  (middle)  and 
June  (bottom)  1990. 
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Fig.  3.3.1.  (cont.)  FINESA  data  recording  uptime  for  July  (top),  August 
(middle)  and  September  (bottom)  1990. 
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3.4 


Event  detection  operation 


NORESS  detections 

The  number  of  detections  (phases)  reported  during  day  091, 1990,  through 
day  273,  1990,  was  43734,  giving  an  average  of  239  detections  per  processed 
day  (183  days  processed). 

Table  3.4.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS., 
Events  automatically  located  by  NORESS 

During  days  091,  1990,  through  273,  1990,  2275  local  and  regional  events 
were  located  by  NORESS,  based  on  automatic  association  of  P-  and  S-type 
arrivals.  This  gives  an  average  12.4  events  per  processed  day  (183  days  pro¬ 
cessed).  60  %  of  these  events  are  within  300  km,  and  86  %  of  these  events  are 
within  1000  km. 

ARC  ESS  detections 

The  number  of  detections  (phases)  reported  during  day  091, 1990,  through 
day  273,  1990,  was  84203,  giving  an  average  of  460  detections  per  processed 
day  (183  days  processed). 

Table  3.4.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 
Events  automatically  located  by  ARCESS 

During  days  091,  1990,  through  273,  1990,  4307  local  and  regional  events 
were  located  by  ARCESS,  based  on  automatic  association  of  P-  and  S-type 
arrivals.  This  gives  an  average  23.5  events  per  processed  day  (178  days  pro¬ 
cessed).  49  %  of  these  events  are  within  300  km,  and  84  %  of  these  events  are 
within  1000  km. 

FINESA  detections 

The  number  of  detections  (phases)  reported  during  day  091,  1990,  through 
day  273,  1990,  was  28235,  giving  an  average  of  392  detections  per  processed 
day  (72  days  processed). 

Table  3.4.3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 
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Events  automatically  located  by  FINESA 

During  days  091, 1990,  through  273, 1990,  1385  local  and  regional  events 
were  located  by  FINESA,  based  on  automatic  association  of  P-  and  S-type 
arrivals.  This  gives  an  average  19.2  events  per  processed  day  (72  days  pro¬ 
cessed).  68  %  of  these  events  are  within  300  km,  and  87  %  of  these  events  are 
within  1000  km. 

J.  Fyen 
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3.5  IMS  operation 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  De¬ 
cember  1989  and  has  been  operated  experimentally  since  1  January  1990  for 
automatic  processing  of  multiple-array  data.  The  current  version  of  IMS  pro¬ 
cesses  data  from  the  two-array  network  consisting  of  NORESS  and  ARCESS. 
Future  upgrades  of  IMS  will  allow  data  from  additional  arrays  and  single  sta¬ 
tions  to  be  incorporated. 

In  general,  our  routine  operation  of  the  IMS  during  the  reporting  period 
has  progressed  well,  and  the  system  has  proved  to  be  very  powerful  and  flex¬ 
ible.  The  well-developed  automatic  processing  combined  with  very  versatile 
interactive  tools  has  kept  the  analyst  workload  at  a  low  level,  and  in  fact  only 
one  analyst  has  been  required  to  handle  the  regular  processing  at  any  time. 

Since  the  IMS  is  still  in  an  initial  stage,  there  have  naturally  been  some 
problems  of  a  technical  nature,  but  these  have  diminished  considerably  as  the 
system  has  matured.  At  the  end  of  the  reporting  period,  the  IMS  system  is 
running  very  smoothly. 

Table  3.5.1  presents  an  overview  of  IMS  event  processor  downtimes.  Table 
3.5.2  gives  a  summary  of  phase  detections  and  processed  regional  events  by 
IMS.  From  top  to  bottom,  the  table  gives  the  total  number  of  detections  by 
the  IMS,  the  detections  that  are  associated  with  regional  events  declared  by 
the  IMS,  the  number  of  detections  that  are  not  associated  with  such  events, 
the  number  of  regional  events  declared  by  the  IMS,  the  number  of  such  events 
rejected  by  the  analyst,  the  number  of  events  accepted  by  the  analyst,  the 
number  of  events  accepted  by  the  analyst  without  any  changes,  and  finally  the 
number  of  events  accepted  after  some  sort  of  modification  by  the  analyst.  This 
last  category  is  divided  into  three  classes:  Events  where  phases  (not  detected 
by  the  IMS)  have  been  added  by  the  analyst,  events  for  which  the  phase 
assignments  by  the  IMS  have  been  changed  or  one  or  more  phase  detections 
have  been  removed,  and  events  for  which  the  changes  by  the  analyst  have 
amounted  to  retiming  the  phases  only. 

From  initial  review  of  the  IMS  bulletin,  it  is  clear  that  the  final  output  is  of 
very  high  quality  from  a  seismological  point  of  view.  More  detailed  evaluation 
will  be  conducted  at  a  later  stage. 
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1990  May  29  14:47:01.000  - 
1990  May  30  13:13:30.000  - 
1990  Jun  8  14:20:32.713  - 
1990  Jun  9  12:55:30.000  - 
1990  Jun  10  11:22:22.752  - 
1990  Jun  15  13:23:21.700  - 
1990  Jun  16  13:24:57.500  - 
1990  Jun  22  20:13:36.000  - 
1990  Jun  26  03:34:17.000  - 
1990  Jun  30  10:20:53.800  - 
1990  Jul  6  07:45:52.000  - 
1990  Jul  9  17:56:50.500  - 
1990  Jul  10  23:29:23.700  - 
1990  Jul  13  13:26:07.600  - 
1990  Jul  16  16:15:52.800  - 
1990  Jul  17  09:17:33.300  - 
1990  Jul  19  21:39:56.200  - 
1990  Jul  22  20:27:00.000  - 
1990  Jul  24  09:55:22.500  - 
1990  Aug  5  05:26:02.840  - 
1990  Aug  13  20:15:28.400  - 
1990  Aug  18  20:38:30.200  - 
1990  Sep  1  05:14:22.800  - 
1990  Sep  10  04:43:55.100  - 
1990  Sep  28  09:40:09.800  - 
1990  Sep  30  10:29:13.300  - 


1990  May  30  10:06:07.000 
1990  May  31  18:30:50.379 
1990  Jun  9  08:00:00.000 
1990  Jun  9  22:30:51.000 
1990  Jun  11  10:00:00.000 
1990  Jun  16  00:00:00.000 
1990  Jun  19  00:00:00.000 
1990  Jun  26  00:00:00.000 
1990  Jun  26  18:27:20.000 
1990  Jul  1  00:00:00.000 
1990  Jul  7  00:00:00.000 
1990  Jul  10  00:00:00.000 
1990  Jul  12  00:00:00.000 
1990  Jul  15  00:00:00.000 
1990  Jul  17  00:00:00.000 
1990  Jul  19  00:00:00.000 
1990  Jul  22  10:56:47.400 
1990  Jul  24  00:24:57.933 
1990  Jul  24  23:35:08.241 
1990  Aug  6  05:18:15.167 
1990  Aug  17  18:00:00.000 
1990  Aug  19  18:00:00.000 
1990  Sep  2  00:00:00.000 
1990  Sep  24  00:00:00.000 
1990  Sep  29  00:00:00.000 
1990  Sep  30  24:00:00.000 


Table  3.5.1.  IMS  event  processor  (pipeline)  downtimes.  Logging  of  time 
intervals  not  processed  by  the  event  processor  (pipeline)  started  the  last  week 
of  May.  The  main  reason  for  these  downtimes  was  slow  performance  of  the 
event  processing,  causing  loss  of  waveform  data  from  the  disk  before  the  events 
were  processed. 
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Apr  90 

May  90 

Jun  90 

Jul  90 

Aug  90 

Sep  90 

total 

Phase  detections 

11843 

13988 

13178 

10989 

19267 

15402 

84667 

-  Associated  phases 

1685 

2521 

1863 

1125 

2361 

1021 

10576 

-  Unassociated  phases 

10158 

11467 

11315 

9864 

16906 

14381 

74091 

Events  declared 

708 

1002 

834 

493 

943 

278 

4258 

-  Rejected  events 

183 

196 

237 

131 

218 

43 

1009 

-  Accepted  events 

525 

806 

597 

361 

725 

235 

3249 

Unchanged  events 

85 

186 

293 

172 

290 

78 

1104 

Modified  events 

440 

620 

304 

189 

435 

157 

2145 

Phases  added 

46 

63 

23 

9 

30 

16 

187 

Phases  changed  or  removed 

281 

402 

270 

180 

405 

141 

1679 

Retiming  only 

113 

155 

11 

0 

0 

0 

279 

Table  3.5.2.  IMS  phase  detections  and  events  summary.  Statist.es  for  “retiming  only” 
have  been  compiled  only  up  to  June  4. 
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4  Improvements  and  Modifications 

4.1  NORSAR 

NORSAR  data  acquisition 

No  modification  has  been  made  to  the  NORSAR  data  acquisition  system. 

RD3  digitizers  and  acquisition  equipment  have  been  acquired  from  Nano¬ 
metrics 

to  test  recording  at  one  subarray.  There  are  two  major  tasks  that  need  to 
be  worked  upon: 

1.  Solving  the  problem  of  getting  enough  DC  current  from  the  central  vault 
to  the  digitizer  at  the  seismometer. 

2.  Including  SDLC  communication  in  the  Nanometrics  system  to  get  a  con¬ 
tinuous  data  stream  from  the  subarray  to  the  NORSAR  data  processing 
center.  (This  is  currently  not  offered  by  the  manufacturer.) 

NORSAR  detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the 
IBM  during  this  reporting  period. 

Detection  statistics  are  given  in  section  2. 

NORSAR  event  processing 

There  are  no  changes  in  the  routine  processing  of  NORSAR  events,  using 
the  IBM  system. 


4.2  NORESS/ARCESS/FINESA/GERESS 

Detection  processing 

The  routine  detection  processing  of  the  arrays  is  running  satisfactorily  on 
each  of  the  arrays’  SUN-3/280  acquisition  systems.  The  same  program  is 
used  for  NORSAR/NORESS/ARCESS/  FINESA/GERESS,  but  with  differ¬ 
ent  recipes.  The  beam  table  for  NORESS/'ARCESS  is  found  in  NORSAR  Sci. 
Rep.  No.  1-89/90.  A  new  beam  table  for  FINESA  is  found  in  Table  4.2.1.  This 
change  of  the  beam  table  is  related  to  changes  in  the  FINESA  array  geome¬ 
try  in  August.  A  new  vertical  sensor  was  installed  at  the  array  center  (AO), 
and  horizontal  sensors  were  added  at  the  existing  A1  site.  Test  processing  of 
GERESS  data  has  been  started,  and  the  initial  beam  table  is  shown  in  Table 
4.2.2. 
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Event  processing.  Phase  estimation 

This  process  does  F-K  and  polarization  analysis  for  each  detection  to  de¬ 
termine  phase  velocity,  azimuth  and  type  of  phase. 

The  processing  using  the  EP  program  is  now  running  on  the  IMS  machines, 
and  detections  are  put  into  the  ORACLE  detection  data  base  for  use  by  IMS. 

Plot  and  epicenter  determination 

A  description  of  single  array  event  processing  is  found  in  NORSAR  Sci. 
Rep.  No. '  -88/89,  and  NORSAR  Sci.  Rep.  No.  2-89/90. 
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BEAM 

Velocity 

Azimuth 

Filter 

C011 

99999.9 

0.0 

0.5  - 

C021 

99999.9 

0.0 

1.0  - 

C031 

99999.9 

0.0 

1.5  - 

C032 

11.0 

30.0 

1.5  - 

C033 

11.0 

90.0 

1.5  - 

C034 

11.0 

150.0 

1.5  - 

C035 

11.0 

210.0 

1.5  - 

C036 

11.0 

270.0 

1.5  - 

C037 

11.0 

330.0 

1.5  - 

C038 

15.0 

80.0 

1.5  - 

C039 

10.0 

20.0 

1.5  - 

C041 

99999.9 

0.0 

2.0  - 

C042 

10.2 

30.0 

2.0  - 

C043 

10.2 

90.0 

2.0  - 

C044 

10.2 

150.0 

2.0  - 

C045 

10.2 

210.0 

2.0  - 

C046 

10.2 

270.0 

2.0  - 

C047 

10.20 

330.0 

2.0  - 

C048 

15.0 

80.0 

2.0  - 

C049 

10.0 

20.0 

2.0  - 

C051 

99999.9 

0.0 

2.5  - 

C052 

8.9 

30.0 

2.5  - 

C053 

8.9 

90.0 

2.5  - 

C054 

8.9 

150.0 

2.5  - 

C055 

8.9 

210.0 

2.5  - 

C056 

8.9 

270.0 

2.5  - 

C057 

8.9 

330.0 

2.5  - 

C058 

15.0 

80.0 

2.5  - 

C061 

99999.9 

0.0 

3.0  - 

C062 

10.5 

30.0 

3.0  - 

C063 

10.5 

90.0 

3.0  - 

C064 

10.5 

150.0 

3.0  - 

C065 

10.5 

210.0 

3.0  - 

C066 

10.5 

270.0 

3.0  - 

C067 

10.5 

330.0 

3.0  - 

C068 

15.0 

80.0 

3.0  - 

C071 

99999.9 

0.0 

3.5  - 

C072 

11.1 

30.0 

3.5  - 

C073 

11.1 

90.0 

3.5  - 

band  Threshold  Configuration 


1.5 

4.4 

7 

C 

3.0 

4.4 

7 

C 

3.5 

4.4 

7 

C 

3.5 

4.4 

7 

C 

3.5 

4.4 

7 

C 

3.5 

4.4 

7 

C 

3.5 

4.4 

7 

C 

3.5 

4.4 

7 

c 

3.5 

4.4 

7 

c 

3.5 

3.9 

7 

c 

3.5 

3.9 

7 

c 

4.0 

4.4 

7 

c 

4.0 

4.4 

7 

c 

4.0 

4.4 

7 

c 

4.0 

4.4 

7 

c 

4.0 

4.4 

7 

c 

4.0 

4.4 

7 

c 

4.0 

4.4 

7 

c 

4.0 

3.9 

7 

c 

4.0 

3.9 

7 

c 

4.5 

4.4 

13 

BC 

4.5 

4.4 

13 

BC 

4.5 

4.4 

13 

BC 

4.5 

4.4 

13 

BC 

4.5 

4.4 

13 

BC 

4.5 

4.4 

13 

BC 

4.5 

4.4 

13 

BC 

4.5 

3.9 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

4.4 

13 

BC 

5.0 

3.9 

13 

BC 

5.5 

4.4 

16 

ABC 

5.5 

4.4 

16 

ABC 

5.5 

4.4 

16 

ABC 
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C074 

11.1 

150.0 

3.5 

- 

5.5 

4.4 

16 

ABC 

C075 

11.1 

210.0 

3.5 

- 

5.5 

4.4 

16 

ABC 

C076 

11.1 

270.0 

3.5 

- 

5.5 

4.4 

16 

ABC 

C077 

11.1 

330.0 

3.5 

- 

5.5 

4.4 

16 

ABC 

C081 

99999.9 

0.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C082 

9.5 

30.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C083 

9.5 

90.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C084 

9.5 

150.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C085 

9.5 

210.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C086 

9.5 

270.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C087 

9.5 

330.0 

4.0 

- 

8.0 

4.4 

16 

ABC 

C091 

92999.9 

0.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C092 

10.5 

30.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C093 

10.5 

90.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C094 

10.5 

150.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C095 

10.5 

210.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C096 

10.5 

270.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C097 

10.5 

330.0 

5.0 

- 

10.0 

4.9 

16 

ABC 

C101 

99999.9 

0.0 

8.0 

- 

16.0 

4.9 

10 

AB 

C102 

9.9 

30.0 

8.0 

- 

16.0 

4.9 

10 

AB 

C103 

9.9 

90.0 

8.0 

- 

16.0 

4.9 

10 

AB 

C104 

9.9 

150.0 

8.0 

- 

16.0 

4.9 

10 

AB 

C105 

9.9 

210.0 

8.0 

- 

16.0 

4.9 

10 

AB 

C106 

9.9 

270.0 

8.0 

- 

16.0 

4.9 

10 

AB 

C107 

9.9 

330.0 

8.0 

- 

16.0 

4.9 

10 

AB 

CI01 

99999.9 

0.0 

0.5 

- 

1.5 

2.9 

7 

C 

CI02 

99999.9 

0.0 

1.0 

- 

2.0 

2.9 

7 

C 

CI03 

99999.9 

0.0 

1.5 

- 

2.5 

2.9 

7 

C 

CI04 

99999.9 

0.0 

2.0 

- 

4.0 

2.5 

7 

C 

CI05 

99999.9 

0.0 

3.5 

- 

5.5 

2.5 

7 

C 

CI06 

99999.9 

0.0 

5.0 

- 

10.0 

2.7 

7 

C 

Table  4.2.1.  FINESA  beam  table.  Valid  from  1990-237:  (25  august  1990), 
when  AO  was  included.  The  table  shows  name  of  beam,  velocity  (km/sec), 
azimuth  (degrees),  filter  band  (Hz),  STA/LTA  threshold,  and  configuration. 
The  configuration  is  described  with  number  of  sensors  and  A,  B  o’.  C-rings. 
The  sensor  at  AO  is  included  in  all  configurations.  CI01  -  CI06  are  incoherent 
beams. 
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BEAM 

Velocity 

Azimuth 

Filter 

G011 

99999.9 

0.0 

0.5  - 

G021 

99999.9 

0.0 

1.0  - 

G031 

99999.9 

0.0 

1.5  - 

G032 

11. 0 

30.0 

1.6  - 

G033 

11.0 

90.0 

1.5  - 

G034 

11.0 

150.0 

1.5  - 

G035 

11.0 

210.0 

1.5  - 

G036 

11.0 

270.0 

1.5  - 

G037 

11.0 

330.0 

1.5  - 

G038 

15.8 

80.0 

1.5  - 

G039 

10.0 

30.0 

1.5  - 

G041 

99999.9 

0.0 

2.0  - 

G042 

10.2 

30.0 

2.0  - 

G043 

10.2 

90.0 

2.0  - 

G044 

10.2 

150.0 

2.0  - 

G045 

10.2 

210.0 

2.0  - 

G046 

10.2 

270.0 

2.0  - 

G047 

10.2 

330.0 

2.0  - 

G048 

15.8 

80.0 

2.0  - 

G049 

10.0 

30.0 

2.0  - 

G05i 

99999.9 

0.0 

2.5  - 

G052 

8.9 

30.0 

2.5  - 

G053 

8.9 

90.0 

2.5  - 

G054 

8.9 

150.0 

2.5  - 

G055 

8.9 

210.0 

2.5  - 

G056 

8.9 

270.0 

2.5  - 

G057 

8.9 

330.0 

2.5  - 

G058 

15.8 

80.0 

2.5  - 

G059 

10.0 

30.0 

2.5  - 

G061 

99999.9 

0.0 

3.0  - 

G062 

10.5 

30.0 

3.0  - 

G063 

10.5 

90.0 

3.0  - 

G064 

10.5 

150.0 

3.0  - 

G065 

10.5 

210.0 

3.0  - 

G066 

10.5 

270.0 

3.0  - 

G067 

10.5 

330.0 

3.0  - 

G068 

15.8 

80.0 

3.0  - 

G069 

10.0 

30.0 

3.0  - 

G071 

99999.9 

0.0 

3.5  - 

G072 

11.1 

30.0 

3.5  - 

G073 

11.1 

90.0 

3.5  - 

band  Threshold  Configuration 


1.5 

4.5 

10 

D 

3.0 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

3.5 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.0 

4.5 

17 

CD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

4.5 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.0 

4.5 

22 

BCD 

5.5 

4.5 

13 

BC 

5.5 

4.5 

13 

BC 

5.5 

4.5 

13 

BC 
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G074 

11.1 

150.0 

3.5  - 

5.5 

4.5 

13 

BC 

G075 

11.1 

210.0 

3.5  - 

5.5 

4.5 

13 

BC 

G076 

11.1 

270.0 

3.5  - 

5.5 

4.5 

13 

BC 

G077 

11.1 

330.0 

3.5  - 

5.5 

4.5 

13 

BC 

G081 

99999.9 

0.0 

4.0  - 

8.0 

4.5 

13 

BC 

G082 

9.5 

30.0 

4.0  - 

8.G 

4.5 

13 

BC 

G083 

9.5 

90.0 

4.0  - 

8.0 

4.5 

13 

BC 

G084 

9.5 

150.0 

4.0  - 

8.0 

4.5 

13 

BC 

G085 

9.5 

210.0 

4.0  - 

8.0 

4.5 

13 

BC 

G086 

9.5 

270.0 

4.0  - 

8.0 

4.5 

13 

BC 

G087 

9.5 

330.0 

4.0  - 

8.0 

4.5 

13 

BC 

G091 

99999.9 

0.0 

5.0  - 

10.0 

5.0 

13 

BC 

G092 

10.5 

30.0 

5.0  - 

10.0 

5.0 

13 

BC 

G093 

10.5 

90.0 

5.0  - 

10.0 

5.0 

13 

BC 

G094 

10.5 

150.0 

5.0  - 

10.0 

5.0 

13 

BC 

G095 

10.5 

210.0 

5.0  - 

10.0 

5.0 

13 

BC 

G096 

10.5 

270.0 

5.0  - 

10.0 

5.0 

13 

BC 

G097 

10.5 

330.0 

5.0  - 

10.0 

5.0 

13 

BC 

G101 

99999.9 

0.0 

8.0  - 

16.0 

5.0 

9 

AB 

G102 

9.9 

30.0 

8.0  - 

16.0 

5.0 

9 

AB 

G103 

9.9 

90.0 

8.0  - 

16.0 

5.0 

9 

AB 

G104 

9.9 

150.0 

8.0  - 

16.0 

5.0 

9 

AB 

G105 

9.9 

210.0 

8.0  - 

16.0 

5.0 

9 

AB 

G106 

9.9 

270.0 

8.0  - 

16.0 

5.0 

9 

AB 

G107 

9.9 

330.0 

8.0  - 

16.0 

5.0 

9 

AB 

GH01 

99999.9 

0.0 

2.0  - 

4.0 

3.0 

8 

AC_ne 

GH02 

99999.9 

0.0 

3.5  - 

5.5 

3.0 

8 

AC_ne 

GH03 

1  '999.9 

0.0 

5.0  - 

10.0 

3.0 

8 

AC_ne 

GH04 

99999.9 

0.0 

8.0  - 

16.0 

3.0 

8 

AC.ne 

GV01 

99999.9 

0.0 

0.5  - 

1.5 

3.0 

10 

D 

GV02 

99999.9 

0.0 

1.0  - 

2.0 

3.0 

8 

C 

GV03 

99999.9 

0.0 

1.5  - 

2.5 

3.0 

8 

C 

GV04 

99999.9 

0.0 

2.0  - 

3.0 

3.0 

8 

C 

GV05 

99999.9 

0.0 

2.0  - 

4.0 

3.0 

8 

C 

GV06 

99999.9 

0.0 

3.5  - 

5.5 

3.0 

8 

C 

Table  4.2.2.  GERESS  beam  table.  Initial  test  beam  set.  The  table  shows 
name  of  beam,  velocity  (km/sec),  azimuth  (degrees),  filter  band  (Hz),  STA/LTA 
threshold,  and  configuration.  The  configuration  is  described  with  number  of 
sensors  and  A,  B,  C  or  D-rings.  Center  AO  is  always  within  the  configura¬ 
tion.  GH01  -  GH04  are  incoherent  beams  using  horizontal  components  of  the 
A2,  Dl,  D4  and  D7  three- component  stations.  GV01  -  GV04  are  incoherent 
beams  using  vertical  components.  This  beam  set  has  been  in  operation  since 
1090-294:03. 
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5  Maintenance  Activities 

5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  in  the  field,  at  the  Maintenance  Center 
(NMC)  Hamar  and  NDPC  activities  related  to  monitoring  and  control  of  the 
NORSAR,  NORESS  and  ARCESS  arrays. 

Activities  involve  preventive  and  corrective  maintenance,  and  in  addition 
installation  and  modification  of  equipment  in  the  field  and  at  the  Maintenance 
Center. 

NORSAR 

Visits  were  made  to  the  arrays  in  connection  v/ith  adjustment  of  channel 
gain  and  offset  SP/LP.  Adjustments  of  LP  seismometer  mass  position  (MP) 
and  free  period  (FP)  were  made.  Cables  have  been  located/pointed  out  and 
also  repaired/spliced.  Subarrays  have  been  visited  in  connection  with  loss  of 
power.  At  the  Well  Head  Vaults  (WHVs)  RA-5  amplifiers  have  been  replaced. 
Communication  lines  have  been  tested  and  modems  replaced. 

NORESS 

The  NORESS  array  has  been  visited  in  connection  with  replacement  of 
fiber  optical  transmitters  on  the  Hub  card  (site  C4).  A  water  leakage  was 
stopped  in  connection  with  C2  WHV.  In  June  an  AC  power  failure  caused  a 
visit  to  the  central  station  (Hub).  In  September  a  new  UPS  (Uninterruptible 
Power  Supply)  system  was  installed. 

ARCESS 

Also  here  activities  have  been  diverse.  An  air-conditioner  in  the  central 
building  ac  the  array  site  has  been  replaced,  and  fiber  optical  transmitters  have 
been  replaced  on  fiber  optical  links.  Fiber  optical  transmitters/receivers  on 
the  Hub  15  card  were  replaced  (site  C4).  In  addition,  transmitters/receivers 
on  a  DHL  70  card  were  replaced  (site  D4).  A  GS-  13  seismometer  was  replaced 
(site  D5). 

Details  are  presented  in  Table  5.1. 


50 


Subarray/ 

Task 

Date 

area 

NORSAR: 

No  visits 

NORESS: 

Replaced  fiber  optical  transmitters  on  Hub  15  card 

5  April 

NMC: 

Tested/modified  new  equipment  in  connection  with 
extension  of  the  FINESA  array 

April 

NDPC: 

Daily  check  of  NORSAR,  NORESS,  ARCESS  and  FINESA 
status 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 

Period.  Adjusted  when  outside  tolerances 

April 

NORSAR: 

02C 

Located  faulty  SLEM  power  supply 

3  May 

01B 

Meausrements  revealed  faulty  cable  SP02,04 

7  May 

01B 

Cable  work  SP02,04 

14  May 

06C 

Channel  gain  and  DC  offset  on  all  SP/LP 
channels  adjusted 

19  May 

01B 

Cable  work  SP02,04  terminated 

22  May 

04C 

SP/LP  checked 

28  May 

02C 

Power  failure  stated;  electricians  from 
the  local  power  plant  had  already  arrived 

30  May 

NORESS: 

Stopped  water  leakage  to  the  C2  Well  Head 

8,9,10 

Vault 

May 

Reinstalled  equipment,  including  seismometers 

10  May 

C4  repaired 

18  May 

NDPC: 

Daily  check  of  NORSAR,  NORESS,  ARCESS  and  FINESA 
status 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 

Period 

May 

Table  5.1  Activities  in  the  field  and  the  NORSAR  Maintenance  Center,  in¬ 
cluding  NDPC  activities  related  to  the  NORSAR,  NORESS,  ARCESS  and 
FINESA  arrays,  1  April  —  30  September  1990. 


Subarray/ 

area 

Task 

Date 

NORSAR: 

02B 

A  modem  was  replaced 

8  June 

AC  failure 

Gain  adjusted  on  all  channels  (LTA) 

DC  offset  adjusted  on  all  channels 

19  June 

“Old”  modem  reinstalled 

22  June 

04C 

Lost  AC  power 

1  June 

06C 

AC  failure,  2  phases  lost 

28  June 

NORESS: 

AC  failure 

28  June 

ARCESS: 

Replaced  air-conditioner 

25,26 

Replaced  optical  transmitter  at  remote  site 

D2 

Replaced  all  fiber  optical  links  (to  remote 
sites) 

Replaced  fiber  optical  transmitters  at  the  Hub 

15  card,  site  B1  and  B3 

June 

NDPC: 

Daily  check  of  NORSAR,  NORESS,  ARCESS  and  FINESA 
status 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 

Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

June 

NORSAR: 

02B 

Visit  in  connection  with  loss  of  power 

12  July 

SLEM  reset 

16  July 

02C 

Replaced  FP/MP  centering  device  EW  seism. 

4  July 

Cable  splicing  carried  out  SP05 

10  July 

04C 

Comm,  work  in  cooperation  with  NTA/Hamar 

3,9  July 

SP02,03  repair 

17  July 

Repair  of  LP  EW  seismometer 

18,19  July 

Gain  and  offset  adjustments  all  channels  SP/LP 

20  July 

06C 

Installed  a  new  RA-5  amplifier  SP02 

12  July 

Cable  work  carried  out  at  site  D5 

16  July 

ARCESS: 

Replaced  fiber  optical  transm. /receiver  on 

Hub  15  card,  site  4 

Reset  “white  box”  communications  unit 

Replaced  transm. /receiver  on  DHL  68  card 
site  AO 

24  July 

Table  5.1  (cont.) 
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Subarray/ 

Task 

Date 

area 

Replaced  DHL  70  card  site  D4  (vertical  channel) 

Replaced  GS13  seismometer  at  site  D5 

NDPC: 

Daily  check  of  NORSAR,  NORESS,  ARCESS  and  FINESA 
status 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 

Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

July 

NORSAR: 

Visit  in  connection  with  power  loss 

27  Aug 

02B 

Cable  work 

27  Aug 

04C 

Repaired  SP  cal.  generator 

27  Aug 

06C 

Cable  pointed  out 

31  Aug 

FINESA: 

NMC  staff  reinstalled  data  acquisition  unit 

20-26 

New  channels  AOZ,  AOZ  low  channel,  A1NS,  AlEW 
installed 

Aug 

NDPC: 

Daily  check  of  NORSAR,  NORESS,  ARCESS  and  FINESA 
status 

Weekly  calibration  of  NORSAR  SP/LP  instruments 
Continuous  measurement  of  Mass  Position  and  Free 

Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Aug 

NORSAR: 

01A 

Replaced  RA-5,  SP02 

Channel  gain  all  SP  channels  adjusted 

MP/FP  NS  and  EW  LP  seismometer  adjusted 

27  Sep 

04  C 

Replaced  RCD  (Remote  Centering  Device)  NS 

LP  seismometer 

11  Sep 

Repaired  improper  connection  in  the  EW  data 
circuit 

06C 

Cable  SP05  pointed  out 

7  Sep 

Cable  work  SP04 

10  Sep 

Splicing  of  cable  SP04,05 

24  Sep 

NORESS: 

Installed  new  UPS  system 

19-20  Sep 

Installation  and  initialization 

26  Sep 

Table  5.1  (cont.) 
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Subarray/ 

Task 

Date 

area 

NDPC: 

Dally  check  of  NORSAR,  NORESS,  ARCESS  and  FINESA 
status 

Weekly  calibration  of  NORSAR  SP/LP  instruments. 
Continuous  measurement  of  Mass  Position  and  Free 

Period 

Adjustment  of  Mass  Position  and  Free  Period  when 
outside  tolerances 

Sep 

Table  5.1  (cont.) 


5.2  Array  status 

As  of  30  September  1990  the  following  NORSAR  channels  deviated  from  tol¬ 
erances. 

01A  01  8  Hz  filter 
02  8  Hz  filter 
04  30  dB  attenuator 


O.A.  Hansen 
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7  Summary  of  Technical  Reports  /  Papers  Pub¬ 
lished 

7.1  Detection  and  yield  estimation  studies 

1.  Yield  estimation  using  Lg  recordings 

Over  the  past  several  years,  extensive  research  has  taken  place  at  NORSAR 
to  develop  and  evaluate  the  RMS  Lg  estimation  technique  for  underground 
nuclear  explosions.  Some  of  the  most  recent  results  have  been  documented  by 
Ringdal  and  Marshall  (1989)  and  Hansen,  Ringdal  and  Richards  (1990). 

These  studies  have  so  far  concentrated  on  the  Shagan  River  (Semipalatinsk) 
test  site,  using  data  from  stations  in  Norway,  Germany,  USSR  and  China  (Fig. 
7.1.1).  By  using  NORSAR  as  a  reference  system,  and  plotting  observed  RMS 
Lg  for  other  f cations  against  NORSAR  values,  it  has  been  found  that  the 
standard  deviations  of  the  residuals  are  consistently  as  low  as  0.03  magnitude 
units,  as  shown  for  selected  stations  in  Fig.  7.1.2.  Furthermore,  comparing 
two  of  the  stations  with  lowest  Lg  detection  threshold  (GAM  and  ARU  in  the 
Soviet  Union),  we  have  found  that  this  consistency  appears  to  span  two  full 
orders  of  event  magnitude,  down  to  approximately  m;,  =  4.0  (Hansen  et  al, 
1990). 

A  possibility  to  compare  the  RMS  Lg  magnitudes  to  published  yields 
for  Semipalatinsk  explosions  has  now  emerged  with  the  recent  publication 
by  Soviet  scientists  quoting  yield  estimates  for  a  number  of  such  explosions 
(Bocharov  et  al,  1989;  see  also  Vergino,  1989a, b).  We  have  made  an  effort  to 
conduct  detailed  analysis  of  available  NORSAR  Lg  data  for  the  explosions  in 
this  data  set.  Since  NORSAR  is  located  at  more  than  4000  km  distance  from 
Semipalatinsk,  only  explosions  of  more  than  10  kt  yield  will  usually  provide 
sufficient  SNR  for  Lg  measurements  at  this  station,  and  even  so,  it  must  be 
noted  that  at  these  lower  yields,  the  estimated  RMS  Lg  is  less  precise  than  at 
yields  exceeding  50  kt. 

Table  7.1.1  gives  the  NORSAR  Lg  magnitudes  obtained  in  these  analyses. 
The  table  comprises  altogether  8  explosions  from  both  the  Shagan  River  (Bal- 
apan),  Degelen  Mountains  and  Konystan  (Murzhik)  subareas.  The  NORSAR 
data  for  the  four  largest  of  these  explosions  have  been  previously  presented  in 
Ringdal  (1989). 

Fig.  7.1.3  (top  part)  shows  a  plot  of  NORSAR  M(Lg)  versus  published 
yields  for  these  8  explosions.  For  comparison,  the  bottom  part  of  Fig.  7.1.3 
shows  world-wide  mb  (UK  values  as  quoted  by  Vergino,  1989a)  versus  yield  for 
the  same  8  events.  The  slope  in  each  plot  has  been  restricted  to  0.75,  and  the 
standard  deviation  in  the  vertical  direction  is  0.043  for  M(Lg)  and  0.081  for  m*, 
versus  log  yield.  We  note  that  this  data  set  is  too  small  to  allow  any  confident 
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estimate  of  the  accuracy  in  estimating  yield  from  M(Lg),  especially  in  view 
of  the  aforementioned  uncertainty  in  M(Lg)  values  for  the  smaller  explosions. 
Nevertheless,  the  figure  would  appear  to  confirm  the  potential  of  RMS  Lg  as 
a  stable  estimator  of  yields  for  fully  coupled  explosions  from  this  region,  and 
indicates  that  even  single  station  M(Lg)  is  better  correlated  to  log(yield)  than 
is  world-wide  mj. 

2.  P-wave  detectability  studies 

Ringdal  (1990)  has  conducted  a  study  of  the  NORESS  array  detection  ca¬ 
pability  for  Semipalatinsk  explosions,  both  in  terms  of  m*,  and  yield  (Bocharov 
et  al,  1989).  In  terms  of  NORESS  mj,  the  50  per  cent  detection  threshold  is 
estimated  at  mj,  =  3.7  ±  0.1.  A  noteworthy  feature,  illustrated  in  Fig.  7.1.4,  is 
the  large  difference  in  NORSAR  mj  bias  (S)  for  the  subregions  Shagan  River 
(S  ss  1.0)  and  Degelen/Konystan  (S  «  0.4).  In  terms  of  world-wide  m^,  the 
estimated  NORESS  detection  threshold,  at  the  50  per  cent  level,  thus  becomes 
m b  —  2.7  and  mj  =  3.3  for  these  subareas,  respectively. 

Fig.  7.1.5  shows  NORESS  m j,  values  (measured  at  the  NORSAR  seismome¬ 
ter  06C02  which  is  co-located  with  the  present  NORESS  center  site)  versus 
log  yield  for  6  low-yield  nuclear  explosions  (the  only  ones  in  Bocharov  et  al 
(1989)  with  available  data  that  did  not  exceed  NORSAR’s  dynamic  range). 
As  before,  we  have  used  a  restricted  slope  of  0.75  in  the  magnitude-log  yield 
relationship,  and  fitted  a  straight  line  with  this  slope  to  the  observed  data. 
The  NORESS  beam  threshold  range  (plus/minus  two  standard  deviations)  is 
shown  as  dotted  horizontal  lines. 

Fig.  7.1.5  indicates  that,  given  similar  coupling  conditions  and  noise  levels 
as  in  our  data  base  of  6  explosions,  Semipalatinsk  explosions  of  yields  at  0.1 
kt  would  be  expected  to  produce  detectable  signals  at  NORESS.  We  note  that 
the  6  reference  explosions  are  all  from  Degelen  or  Konystan.  In  view  of  the 
previous  discussions,  the  detectability  for  Shagan  River  explosions  would  be 
expected  to  be  even  better. 

In  Fig.  7.1.6,  we  plot  NORESS,  ARCESS,  FINESA  and  GERESS  record¬ 
ings  of  an  mi  =  5.7  nuclear  explosion  at  Novaya  Zemlya.  The  figure  shows  an 
unfiltered,  single  sensor  trace  for  each  array,  and  the  signal-to-noise  rato  on 
the  best  array  beam  (STA/LTA  as  computed  by  the  online  detector)  is  given 
in  the  figure  caption.  ARCESS  has  a  strong  S-phase  detection,  and  there  is 
some  S-wave  energy  visible  also  for  the  three  other  arrays.  The  Lg  phase  is  not 
visible  on  these  unfiltered  recordings,  but  in  the  best  filt.pr  band,  the  energy 
in  the  Lg  window  still  exceeds  significantly  the  pre-P  noise. 

From  Fig.  7.1.6  it  is  seen  that  the  P-wave  detection  capability  for  Novaya 
Zemlya  explosions  is  excellent  for  NORESS,  FINESA  and  ARCESS,  in  par¬ 
ticular  the  latter.  Thus,  ARCESS  SNR  for  the  event  shown  was  8383  on  the 
rrray  beam  in  the  filter  band  3-5  Hz.  This  is  more  than  3  orders  of  magni- 
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tude  above  the  operational  threshold,  and  we  note  that  similar  SNR  has  been 
observed  for  the  other  Novaya  Zemlya  explosions  of  similar  size  recorded  at 
ARCESS  (only  4  such  events  have  been  recorded  since  the  array  came  into 
operation).  While  this  is  too  limited  a  data  base  to  give  a  reliable  threshold 
estimate,  a  straight  extrapolation  would  indicate  that  the  ARCESS  threshold 
at  this  site  is  well  below  mj  =  3.0.  NORESS  and  FINESA  also  appear  to  have 
thresholds  close  to  mi  =  3.0,  whereas  GERESS  has  a  lower  detectability  for 
this  test  site. 

3.  Conclusions  and  Recommendations 

Our  studies  confirm  that  Lg  magnitude  estimates  of  Semipalatinsk  explo¬ 
sions  are  remarkably  consistent  between  stations  widely  distributed  in  epicen- 
tral  distance  and  azimuth.  It  thus  appears  that  a  single  station  with  good 
signal- to- noise  ratio  can  provide  M(Lg)  measurements  with  an  accuracy  (one 
standard  deviation)  of  about  0.03  magnitude  unit.  It  is  noteworthy  that  this 
accuracy  is  consistently  obtained  for  a  variety  of  stations  at  very  different 
azimuths  and  distances,  even  though  the  basic  parameters  remain  exactly  as 
originally  proposed  by  Ringdal  for  NORSAR  recordings  (0. 6-3.0  Hz  bandpass 
filter,  RMS  window  length  of  2  minutes,  centered  at  a  time  corresponding  to 
a  group  velocity  of  3.5  km/s).  Moreover,  the  Lg  phase  shows  considerable 
promise  for  use  in  yield  determination,  although  more  data  will  be  needed 
before  the  accuracy  of  Lg-estimated  yields  can  be  firmly  established. 

The  excellent  detection  capability  of  the  regional  arrays  in  northern  Eu¬ 
rope  have  been  confirmed  by  case  studies,  comprising  a  detailed  evaluation  of 
the  NORESS  capability  for  Semipalatinsk  explosions,  and  preliminary  obser¬ 
vations  of  ARCESS  detection  of  Novaya  Zemlya  explosions. 

Further  work  will  be  directed  toward  additional  expansion  of  these  stud¬ 
ies,  in  particular  analyzing  new  data  as  they  become  available,  and  pursuing 
our  analysis  of  Novaya  Zemlya  recordings  at  regional  arrays  as  well  as  other 
available  stations. 


F.  Ringdal 
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Date 

Region 

mt 

NORSAR 

M(Lg) 

Yield 

(kt) 

11/30/69 

Shagan  (TZ) 

6.048 

6.043 

125 

04/25/71 

Degelen 

6.076 

5.862 

90 

06/06/71 

Konystan 

5.526 

5.44 

16 

10/09/71 

Konystan 

5.371 

5.25*) 

12 

10/21/71 

Konystan 

5.580 

5.54 

23 

02/10/72 

Shagan  (NE) 

5.370 

5.37 

16 

11/02/72 

Shagan  (SW) 

6.224 

6.118 

165 

12/10/72 

Shagan  (NE) 

5.996 

6.095**) 

140 

*)  Low  precision  due  to  low  SNR 
**)  Adjusted  for  interfering  explosion 


Table  7.1.1.  Parameters  of  eight  explosions  analyzed  in  this  study.  Yield 
estimates  are  from  Bocharov  et  al  (1989),  and  the  m^  values  are  the  “UK  m 
listed  in  Vergino  (1989a). 
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Fig.  7.1.1.  Map  sh< 
(triangles).  The  two  m 
are  marked  as  filled  ci 


Station  distribution 


ving  the  location  of  stations  analyzed  in  this  study 
in  USSR  test  sites  (Novaya  Zemlyaand  Semipaiatinsk) 


RMS  Lg  comparison 


a) 

Z  Component  RMS  Lg  Comparison 

S=0.92  1=  0.09  S.D.=  0.035  N=  5 


e) 

Z  Component  RMS  Lg  Comparison 

S=1.04  1=  0.71  S.D.=  0.024  N=  12 


b) 

Z  Component  RMS  Lg  Comparison 

S=0.96  1=  1.20  S.D.=  0.025  N=  5 


d) 

Z  Component  RMS  Lg  Comparison 

S=1.48  1=  -2.24  S.D.=  0.023  N=  9 


Fig.  7.1.2.  Examples  of  observed  correlations  between  RMS  Lg  at  NORSAR 
and  selected  stations  in  the  USSR  and  China.  Note  the  excellent  consistency 
in  all  of  the  plots. 
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MAGNITUDE  -  YIELD 

SEMIPALA  TINSK  EXPLOSIONS 
S=  0.75  1=  4.470  D=  0.043  N=  8 


NORSAR  M(Lg) 


I  I  1  I  I  I  »  * 

T 1.0  1.2  1.4  1.6  1.8  2.0  2.2  2.4 

Log  Yield  (Bocharov  et  al) 


MAGNITUDE  -  YIELD 

SEMIPALATINSK  EXPLOSIONS 
0.75  1 =  4.530  D=  0.081  N= 


World-wide  mb 

/' 

/' 

s' 

t 

Log  Yield  (Bocharov  et  al) 


Fig.  7.1.3.  Magnitude-yield  relationships  for  8  nuclear  explosions  at  Semi- 
palatinsk,  with  yields  provided  by  Bocharov  et  al  (1989):.  a)  NORSAR  M(Lg) 
versus  log(yield)  and  b)  world-wide  mi  versus  log(yield).  The  slopes  of  the 
straight  lines  have  been  restricted  to  0.75.  Note  that  the  NORSAR  values 
show  significantly  less  scatter. 


Fig.  7.1.4.  Comparison  of  NORESS  and  world-wide  m &  (as  calculated  at 
Blacknest,  United  Kingdom)  for  Semipalatinsk  explosions.  Note  the  difference 
in  average  m6  bias  between  events  from  Shagan  River  (1.0  m&  units)  and 

Dogolon  Mountains  (0.4  mt  units).  The  straight  lines  on  the  plot  have  a 
restricted  slope  of  1.00. 


Fig.  7.1.5.  Observed  mj,  versus  yield  (Bocharov  et  al,  1989)  for  six  Semi- 
palatinsk  nuclear  explosions  listed  in  Table  7.1.1.  The  straight  line  has  been 
fitted  using  a  restricted  slope  of  0.75.  The  mj  values  are  based  on  NORSAR 
seismometer  06C02,  located  at  the  present  NORESS  center  site.  The  esti¬ 
mated  range  of  the  NORESS  m^  detection  threshold  is  indicated  (see  text  for 
details). 
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Novaya  Zemtya  nuclear  explosion  24  October  1990  -  Regional  array  recordings 


4132.0 


Fig.  7.1.6.  Selected  single  seismometer  SPZ  recordings  from  the  four  regional 
arrays  in  northern  Europe  for  the  m j,  =  5.7  nuclear  explosion  at  Novaya  Zemlya 
on  24  October  1990.  Epicentral  distances  and  maximum  P-phase  signal-to- 
noise  ratios  on  the  array  beam  are  as  follows: 

ARCESS  1110  km,  SNR  =  8383;  FINESA  1780  km,  SNR  =  2189 
NORESS  2270  km,  SNR  =  2478;  GERESS  3380  km,  SNR  =  105 
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7.2  Initial  results  from  real-time  processing  of  GERESS  ar¬ 
ray  data 

Introduction 

The  GERESS  array  (Harjes,  1990)  represents  the  most  recent  addition  to 
the  network  of  regional  seismic  arrays  in  Europe.  GERESS,  which  is  located 
in  the  Bavarian  Forest  area  in  southern  Germany,  has  a  geometry  similar  to 
that  of  NORESS  in  southern  Norway  and  ARCESS  in  northern  Norway,  but 
with  a  slightly  enlarged  instrument  spacing  (Fig.  7.2.1).  Thus  the  diameter  of 
GERESS  is  4  km,  compared  to  3  km  for  the  two  arrays  in  Nor  way. 

Real-time  data  from  GERESS  are  transmitted  continuously  both  to  Bochum, 
Germany  (by  land  line),  and  to  the  NORSAR  Data  Processing  Center  (NDPC) 
at  Kjeller  (by  satellite).  At  NDPC,  the  data  are  currently  subjected  to  continu¬ 
ous  real-time  detection  processing,  using  the  RONAPP  algorithm  (Mykkeltveit 
and  Bungum,  1984).  In  the  near  future,  these  data  will  also  be  incorporated 
into  the  Intelligent  Monitoring  System  (IMS)  (Bache  et  al,  1990). 

GERESS  real-time  processing 

The  initial  processing  of  GERESS  data  is  conducted  at  NDPC  using  a 
beam  set  essentially  identical  to  that  being  used  for  NORESS  and  ARCESS 
data  (see  Subsection  3.4).  Although  several  modifications  will  undoubtedly 
be  desirable  after  more  experience  is  gained,  it  is  still  of  interest  to  observe 
how  a  “generic”  parameter  set  performs  when  applied  to  a  regional  array  in  a 
new  geological  environment.  Below,  we  give  some  initial  examples  from  this 
processing. 

Table  7.2.1  shows  excerpts  from  the  automatic  GERESS  phase  detection 
list  for  day  293  (23  October  1990).  For  each  phase,  the  detection  time,  beam 
number,  SNR,  phase  velocity  and  azimuth  are  listed,  together  with  a  number 
of  attributes  used  in  further  characterization  of  the  detected  phase. 

In  Fig.  7.2.2  we  show  an  automatically  generated  event  plot  (Fyen,  1989) 
for  one  of  the  located  events  during  the  time  period  covered  by  the  detection 
list  in  Table  7.2.1.  This  event  is  located  in  the  Lubin  mining  district  in  Poland, 
at  a  distance  of  344  km  from  GERESS,  and  is  probably  one  of  a  large  number  of 
rockbursts  known  to  occur  in  that  region  (Gibowics,  1987).  We  note  that  this 
event  was  too  small  to  be  detected  by  either  NORESS,  ARCESS  or  FINESA. 

Fig.  7.2.3  shows  GERESS  automatic  output  plot  from  the  P-wave  detection 
of  the  Novaya  Zemlya  explosion  24  October  1990  (A  =  30.4  degrees).  The 
phase  velocity  is  indicative  of  a  teleseismic  signal,  and  the  estimated  azimuth 
is  consistent  with  the  actual  epicenter.  GERESS  SNR  is  105  on  the  array 
beam,  which  is  naturally  lower  than  for  the  three  other  arrays,  all  of  which 
are  within  regional  distances  from  Novaya  Zemlya. 
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An  interesting  feature  of  the  GERESS  recordings  of  this  event  is  a  very 
sharp  PcP  detection  at  2  min  59.5  seconds  after  P  (Fig.  7.2.4).  The  time 
differential  between  P  and  PcP  gives  a  very  sharp  constraint  on  the  epicentral 
distance,  in  particular  when  the  phase  onset  times  can  be  read  as  precisely  as  in 
this  case.  Thus,  given  appropriate  calibration,  very  useful  location  information 
could  be  derived. 

Conclusions 

Our  initial  experience  shows  that  real-time  processing  of  GERESS  data  will 
form  a  very  valuable  contribution  to  the  regional  array  monitoring  research 
network.  So  far,  it  appears  that  many  new  events  in  southern  Europe  within 
500  km  distance  from  GERESS  will  be  added,  and  GERESS  data  will  provide 
further  constraints  on  the  events  detected  and  located  by  more  than  one  array. 
The  teleseismic  potential  of  GERESS  is  also  quite  promising. 


J.  Fyen 
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EPX  190  1990-293:1 
FKX  376  Phase  PN  at 

6.27.35.996  Lot  51.40  Lon  16.37  Azi  32.56  Ost  344.31  A  1282.0  ,  1282.0  f  1.47  ML  3.91 

1 990-293:16.28.25.875  vel  6.80  ozi  33.30  row  0.190  a  3  o  1282.0  f  4.82  fband  3.50  5.50 
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Fig.  7.2.2.  Automatically  generated  event  plot  for  GERESS.  The  traces 
shown  arc  a  selection  of  single  instruments  filtered  in  the  band  3.5-5.5  Hz.  The 
bottom  two  traces  are  the  P-  and  S-beams  using  the  velocities  and  azimuths 
determined  by  broadband  F-K  analysis.  (The  S-beam,  which  is  the  lowermost 
trace,  is  filtered  in  the  band  1.0-2.0  Hz).  The  time  window  is  6  minutes. 


72 


FKX  1388  beom  G032  At  1990-297:15.04.13.725  snr  105.20 

Phose  P  vel  19.70  ozi  42.30  row  0.510  q  1  o  3792.1  f  2.28.fbond  1.50  3.50 
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Fig.  7.2.3.  Automatically  generated  teleseismic  phase  plot  for  GERESS,  cor¬ 
responding  to  the  P-wave  detection  of  the  Novaya  Zemlya  explosion  on  24 
October  1990.  The  traces  shown  are  a  selection  of  single  instruments  filtered 
in  the  band  1.5-3.5  Hz.  The  bottom  two  traces  are  the  P-beam  filtered  and 
unfiltered,  using  the  velocity  and  azimuth  determined  by  broadband  FK  anal¬ 
ysis.  The  time  window  is  60  seconds. 


FKX  1390  beam  G031  At  1990-297:15.07.13.225  snr  16.0 
Phase  P  vel  41.50  ozi  17,70  rpw  0.560  q  1  a  3792.1  f  1.97.fbond  1,50  3.50 
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Fig.  7.2.4.  Similar  to  Fig.  7.2.3,  but  showing  the  automatic  plot  for  the  PcP 
phase  of  the  Novaya  Zemlya  explosion.  Note  the  clear  onset  and  the  significant 
noise  suppression  on  the  filtered  beam  (second  trace  from  bottom).  Also  note 
that  the  high  estimated  phase  velocity  (41.50  km/s)  is  a  clear  indication  that 
this  is  indeed  a  core  phase. 
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7.3  Generalized  beamforming  using  a  network  ot  four  re¬ 
gional  arrays 

The  generalized  beamforming  method  (GBF)  for  automatic  phase  association 
and  event  location  (Ringdal  and  Kvaerna,  1989)  has  now  been  implemented  for 
processing  of  data  from  the  four  regional  arrays,  ARCESS,  FINESA,  GERESS 
and  NORESS.  This  method  works  from  a  list  of  phase  detections  for  each 
array.  For  a  large  set  of  hypothetical  event  locations,  GBF  searches  for  a 
pattern  of  detections  that  fit  the  theoretically  expected  phase  arrivals  from 
the  hypothetical  locations.  When  a  group  of  matching  detections  is  found, 
the  event  location  having  the  best  fit  to  the  data  is  chosen  as  the  most  likely 
epicenter. 

In  the  following  we  present  a  brief  report  on  our  initial  experience  from 
applying  GFB  experimentally  to  the  four-array  network  during  a  continuous 
five-day  period. 

Data  analysis 

A  total  of  600  beam  grid  points  with  an  average  separation  of  150  km 
is  used  to  span  the  geographical  area  of  interest  (see  Fig.  7.3.1).  This  grid 
represents  the  initial  (coarse)  set  of  test  locations,  and  whenever  a  match  is 
found,  a  b'ampacking  algorithm  with  a  much  denser  grid  is  applied  to  refine 
the  locations. 

The  GBF  output  from  a  typical  12-hour  period  is  given  in  Table  7.3.1. 
For  each  associated  event,  the  table  lists  the  individual  phase  detections  that 
either  have  defined  the  event  or  been  associated  as  coda  phases  (p  or  s).  Phase 
attributes  (SNR,  apparent  velocity,  dominant  frequency,  etc.)  are  also  listed 
for  each  entry. 

From  Table  7.3.1  we  find  that  the  majority  of  the  events  are  observed  by 
one  array  only,  whereas  larger  events  have  phase  observations  at  two  or  more 
arrays.  This  is  consistent  with  earlier  studies  of  two-array  and  three-array 
networks. 

It  is  interesting  to  notice  that  the  incorporation  of  data  from  GERESS  en¬ 
ables  us  to  obtain  more  reliable  event  locations  for  events  in  northern  Europe, 
as  seen  for  the  event  in  Poland  with  origin  time  292:21.57.40  and  epicenter 
location  51.70°N,  15.68:E.  This  event  is  presumed  to  be  a  rockburst  in  the 
Lubin  mining  area,  and  the  location  estimate  is  very  close  to  that  site. 

In  order  to  evaluate  the  performance  of  the  method,  we  have  analyzed 
GBF  results  from  five  consecutive  working  days  (10/22/90  -  10/26/90).  The 
results  can  be  summarized  as  follows: 
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1.  Several  hundred  regional  events  were  associated,  where  the  large  major¬ 
ity  were  detected  by  one  array  only. 

2.  A  total  of  37  events  had  four  or  more  associated  phase  detections,  and 
the  corresponding  locations  are  shown  in  Fig.  7.3.2.  After  inspection 
we  found  that  all  but  one  of  these  events  were  properly  associated  and 
located.  From  Fig.  7.3.2  we  see  that  the  large  majority  of  these  events 
are  located  in  the  active  mining  areas  of  western  USSR.  We  also  notice 
a  precise  location  of  a  nuclear  explosion  at  Novaya  Zemlya. 

3.  The  occurrence  of  incorrect  phase  associations  increases  when  we  allow 
fewer  than  four  phases  to  define  an  event. 


Future  work 

From  inspecting  phase  attributes,  the  analyst  can  in  many  cases  identify 
erroneous  phase  associations  based  on  inconsistencies  among  parameters  like 
epicentral  distance,  apparent  velocity  and  dominant  frequency.  Our  next  step 
will  be  to  compile  statistics  on  such  attributes  and  make  this  kind  of  informa¬ 
tion  available  to  the  GBF  program  for  use  in  the  automatic  computations. 

We  have  also  found  that  there  is  a  need  to  refine  the  travel  time  models  and 
also  to  improve  the  estimation  of  onset  time  for  the  different  phase  detections. 
Thereby,  the  automatic  event  locations  will  become  even  more  precise. 

Dynamic  consistency  checks  using  thresholds  computed  by  the  threshold 
monitoring  method  are  also  expected  to  help  in  identifying  and  removing  un¬ 
likely  phase  associations.  This  will  be  particularly  important  at  the  time  when 
detection  data  from  the  two  three-component  systems  in  Poland  will  be  incor¬ 
porated  into  GBF. 

T.  Kvaerna 
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328.6 

Pn 

12.05.09.5 

-0.6 

327.6 

-1.0 

7.4 

13.2 

4020.5 

9  03 

3 

181 

GER 

217  0 

328.6 

P 

12.05.17.1 

335.9 

7.3 

7  9 

4.5 

831.1 

2.69 

2 

182 

GER 

217.0 

328.6 

s 

12.05.30.5 

330.1 

1.5 

3  3 

12.3 

2018.7 

1.25 

1 

-2 

183 

GER 

217.0 

328.6 

Lg 

12.05.37.0 

0.2 

328.5 

-0.1 

3.9 

3.8 

1857.8 

2  28 

2 

184 

Origin  i 

Live 

Lat  Lon 

Azres  Tiaxei. 

Wres 

Nphase 

1990 

-292:12 

.06.12.0 

59 . 50  24.04 

3. 

52  0.32 

1.20 

3 

Sta 

Dist 

hz 

Ph 

Tiae 

Tres 

Azia 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq 

Pol 

Arid 

FIN 

244.1 

208.2 

Pn 

12.06.48.9 

-0  6 

203  1 

-5.1 

7.1 

28  7 

754 .6 

S .  48 

3 

-1 

115939 

FIN 

244.1 

208.2 

Lg 

12.07.20.2 

-0  1 

206.2 

•2.0 

4.1 

17.0 

1227.  7 

3  23 

3 

-2 

115941 

NRS 

708.0 

95.7 

Lg 

12.09.29.8 

-0.2 

99.1 

3,4 

4 .  j 

3.6 

1032.6 

1  50 

1 

-3 

115946 

NRS 

708.0 

95.7 

s 

12.09.32.7 

96.0 

0  3 

3  7 

2.7 

834.6 

2  12 

2 

115947 

Origin 

ti*e 

Lat  Lon 

Azres  Tiares 

Wres 

Nphase 

1990 

1-292  12 

08.49.0 

59  20  27  52 

2. 

89  0.35 

1  07 

2 

Sta 

Dist 

hz 

Ph 

Time 

Tres 

Aria 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq 

Pol 

Arid 

FIN 

?62.6 

161.7 

Pn 

12  09  28  8 

0.0 

157  8 

-3  9 

7  3 

6.1 

306.5 

9  90 

3 

115948 

FIN 

262.6 

161.7 

Lg 

12.10  03.1 

0  7 

159  8 

-1  9 

3.8 

4.7 

708.1 

7  40 

3 

115949 

Origin 

tine 

Lat  Lon 

Azres  Tiares 

Wres 

Nphase 

1990 

i-292  12.28.03  0 

67.60  32  68 

2 

05  l.ll 

1  63 

3 

Sta 

Dist 

hz 

Ph 

Tiae 

Tres 

Azia 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq  Pol 

Arid 

ARC 

363.9 

123.0 

Pn 

12.28.56.4 

1  1 

121.8 

-1  2 

8.0 

90  0 

3292.0 

4  73 

2 

115970 

ARC 

3C3.9 

123.0 

P 

12  29.04  5 

125  5 

2  5 

7  1 

9.6 

1945.9 

6  38 

3 

115971 

ARC 

363  9 

123  r 

Sn 

12.29.35  2 

-1  8 

118.3 

-4.7 

3.7 

2.5 

1338  2 

4.72 

3 

115973 

ARC 

363.9 

123.0 

Lg 

12.29.44  3 

-0.4 

122.8 

-0.2 

3.2 

3.8 

2015.0 

3  25 

2 

-2 

115974 

ARC 

363  9 

123.0 

s 

12  29  47  5 

125.3 

2.3 

4  3 

3.6 

3113  7 

4  44 

2 

115975 

Origin  tiae 

Lat  Lon 

Azres  Tiares 

Wres 

Nphase 

1990-292:12 

.28  38.0 

50.20  14.30 

1. 

97  3.22 

3  72 

3 

Sta 

Dist 

hz 

Ph 

Tiae 

Tres 

Azia 

Ares 

vel 

Snr 

Aap 

Freq 

Fkq 

Pol 

Arid 

GER 

156.8 

15.7 

Pg 

12  29  04  3 

1.0 

16  6 

0  9 

7.5 

4  6 

2374.7 

9  98 

4 

187 

GFR 

156.8 

15.7 

Rg 

12.29  22.5 

-8  7 

12  3 

-3.4 

3  1 

22  7 

650  6 

1.01 

1 

3 

188 

NRS 

118$  2 

170.4 

Sn 

12.33.08.1 

0  1 

172.0 

1.6 

4.3 

2.5 

409.1 

4.49 

3 

115982 

NRS 

1186.2 

170  4 

s 

12.33  12  3 

173.7 

3.3 

3  9 

11.9 

1522  1 

3  95 

2 

-1 

115985 

Origin 

tine 

Lat  Lon 

Azres  Tiares  Wres 

Nphase 

1990-29 2  12 

.30.18.0 

50  80  12.59 

8. 

,90  4.50 

6  72 

3 

Sta 

Dist 

hz 

Ph 

Tiae 

Tres 

Azia 

Ares 

vel 

Snr 

Aap 

Freq 

Fkq 

Pol 

Arid 

GER 

231.8 

340.2 

Pn 

12.30  55  6 

1.6 

32S.8 

14  4 

6  3 

7  7 

3574  5 

8  32 

3 

190 

GER 

231.8 

340.2 

Lg 

12.31.11.1 

-11  8 

330.2 

10  0 

3  0 

8  6 

1874  4 

1  60 

1 

-3 

192 

NRS 

1108  7 

176.2 

Pn 

12.31.40  9 

-0.1 

173  8 

-2.4 

6.7 

5.3 

280  3 

4  13 

3 

115978 

Table  7.3.1.  GBF  output  from  a  typical  12-hour  period.  For  each  associated  ! 

event,  the  table  lists  the  individual  phase  detections  that  either  have  defined  j 

the  event  (Pn,  Sn,  Lg,  Rg)  or  been  associated  as  coda  phases  (p  or  s).  Phase  ? 

attributes  (SNR,  apparent  velocity,  dominant  frequency,  etc.)  are  also  listed 
for  each  entry.  (Page  1  of  3)  ! 
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Origin  time  Lat 

1990-292:12.31  08.0  56.20 


Lon  Azres  Timres  Wrea  Nphase 
26.39  11.07  3.89  6.66  4 


Sta 

Diat 

AX 

Ph 

Time 

Tree 

Azim 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq  Pol 

Arid 

FIN 

584.9 

178.1 

Pn 

12.32.28.5 

1.2 

155.4 

-22.7 

9.4 

4.5 

77.3 

4.56 

3 

115980 

FIN 

584.9 

178.1 

P 

12.32.56  1 

185.0 

6.9 

7  1 

32.1 

1286.9 

3.91 

1 

115984 

FIN 

584.9 

178.1 

sn 

12.33  24.0 

-6.4 

180.3 

2.2 

4.5 

8.7 

767.9 

2.38 

1  -2 

115986 

NRS 

1000.7 

113.8 

Lg 

12.35.55.0 

7.1 

99.4 

-14  4 

3.4 

2.7 

660.2 

2.35 

2 

115991 

GER 

1185.3 

41.6 

Pn 

12.33.39.4 

-0.9 

46.6 

5.0 

6.2 

9.1 

143  9 

'  4.49 

4 

193 

Origin  time 

Lat  Lon 

Axres  Timres  Wrea 

Nphaae 

1990 

•-292:12 

o 

in 

H 

67  60  32.47 

2. 

34  1. 

17  1.75 

4 

Sta 

Diat 

Ax 

Ph 

Time 

Trea 

Azim 

Area 

vel 

Snr 

Amp 

Freq 

Fkq  Pol 

Arid 

ARC 

390.1 

119.8 

Pn 

12.32.40.0 

-0.5 

122.5 

2.7 

7.6 

202.2 

18647.9 

6.78 

1 

115977 

ARC 

390.1 

119.8 

P 

12.32.46.3 

121.6 

1.8 

8.5 

4.2 

1185.3 

3.09 

3 

115979 

ARC 

390.1 

119.8 

P 

12.32.50.5 

117 . 8 

-2.0 

7.4 

2.8 

1568.1 

8.43 

4 

115983 

ARC 

390.1 

119.8 

Sn 

12.33.21.4 

-3.2 

124.8 

5.0 

4  5 

6.9 

5042.9 

2.92 

2  -1 

115987 

ARC 

390.1 

119.8 

Lg 

12  33.34.6 

0.4 

118.2 

-1.6 

3.2 

7.8 

11233.4 

2.86 

2 

115990 

ARC 

390.1 

119.8 

a 

12.33.39.0 

120.4 

0.6 

4.8 

5.4 

10591.9 

3.67 

2 

115994 

NRS 

1302.7 

44.8 

Pn 

12.34.32.1 

0.5 

44.7 

-0.1 

8.4 

4.0 

265.2 

3.35 

3 

115988 

Origin  time  Lat  Lon  Azres  T Imres  Nrea  Nphase 


1990- 

292:12. 

35.44.0 

61.00  29.03 

4. 

27  2.23 

3.29 

8 

Sta 

Diat 

Az 

Ph 

Time 

Trea 

Azim 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq  Pol 

Arid 

FIN 

166.4 

106.0 

Pn 

12.36.09.8 

-1.1 

117.9 

11  9 

6.6 

29.3 

3083  2 

6.99 

3 

-1 

115989 

FIN 

166.4 

106.0 

Lg 

12.36.29.2 

-1.3 

102.0 

-4.0 

4.0 

25.2 

1883.8 

1.85 

1 

-2 

115992 

FIN 

166.4 

106.0 

Rg 

12  36.37  7 

-2.7 

111  3 

5.3 

3.4 

10.9 

4734.5 

1.24 

1 

115993 

NRS 

948.8 

80.6 

Pn 

12.37.46.5 

-1.0 

83.3 

2  7 

11.4 

8.9 

730.2 

6  35 

2 

115995 

NRS 

948.8 

80.6 

Lg 

12.40.07.4 

-2  0 

80.5 

-0  1 

4.6 

2.6 

835.6 

1  36 

2 

115997 

NRS 

948.8 

80.6 

s 

12.40.12  6 

77.5 

-3.1 

4.1 

3.9 

1353  3 

1.90 

1 

116001 

ARC 

966.4 

168.6 

Pn 

12.37  51.7 

2.0 

169  2 

0.6 

9  1 

9.9 

466.4 

5.42 

2 

115998 

ARC 

966  4 

168.6 

Sn 

12.39.27.9 

1  2 

162.5 

-6  1 

4  7 

3  0 

402.9 

4  30 

2 

116000 

ARC 

966.4 

168.6 

Lg 

12.40.07.8 

-6.5 

172.0 

3.4 

3  5 

2.8 

695  1 

2.41 

2 

-3 

116002 

Origin  time  Lat  Lon  Azres  T lares  Kres  Nphaae 

1990-292  13.18  OS.O  59.50  27.85  0  83  0.19  0.39  2 


Sta 

Diat 

Az 

Ph 

Time 

Trea 

Azim 

Area 

Vel 

Snr 

Amp 

Freq 

Fkq 

Pol 

Arid 

FIN 

237.7 

155.0 

Pn 

13.18.41.7 

0  0 

154  4 

-0  6 

6.8 

7  1 

204.3 

8  91 

3 

116032 

FIN 

237  7 

155.0 

Sn 

13.19.12.4 

0.3 

153.9 

-1.1 

3.8 

5.0 

563  3 

7  06 

3 

-2 

116033 

origin 

time 

Lat  Lon 

Azres  Timres  Wres 

Nphaae 

1990-292:13.36.54.0 

59  50  24  83 

7. 

98  0. 

66  2.6S 

3 

Sta 

Diat 

Az 

Ph 

Time 

Tres 

Azim 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq 

Pol 

Arid 

FIN 

227.4 

198.1 

Pn 

13.37.30.0 

0.6 

191  0 

-7.1 

7  8 

16.5 

392.4 

3  37 

1 

116045 

FIN 

227.4 

198  1 

Lg 

13.37.56  3 

-1  4 

193.0 

-5.1 

4.0 

23.2 

1195  6 

2.83 

1 

-2 

116046 

NRS 

750.8 

94.7 

Lg 

13  40  24.1 

0  0 

83  1 

-11.6 

4  1 

2.5 

533.1 

2  58 

2 

llt-050 

Origin  time  Lat  Lon  Azres  Timres  Wres  Nphaae 

1990-292:13.49.09  0  60  40  29.68  i  23  0.47  0.7/,  2 


Sta 

Dist 

Az 

Ph 

Time 

Tres 

Azim 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq 

Pol 

Arid 

FIN 

227.6 

119.2 

Pn 

13  49.44  1 

-0.4 

118.6 

-0.6 

7.1 

6.6 

460.8 

9.0  > 

3 

1 

116071 

FIN 

227  6 

119  2 

a 

13  50.06.4 

121  2 

2  0 

4.1 

5.8 

769  0 

5.57 

3 

-3 

116073 

FIN 

227.6 

119  2 

Lg 

13.50.12  1 

-0  6 

117.3 

-1.9 

3.1 

3.9 

444.4 

1.14 

1 

-3 

116074 

Origin 

time 

Lat  Lon 

Azres  Timres  Wres 

Nphase 

1990 

-292  14 

.24.30.0 

57.40  6  13 

10. 

65  1. 

13  3  79 

3 

Sta 

Diat 

Az 

Pb 

Time 

Tres 

Azim 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq 

Pol 

Arid 

NRS 

484  2 

222.2 

Pn 

14.25.36  3 

-0.8 

211.3 

-10.9 

9.0 

10.8 

635.0 

8  20 

1 

1 

116090 

NRS 

484  2 

222  2 

Lg 

14  26.44  9 

-0.5 

224.7 

2.5 

4  0 

3.6 

909  2 

2  13 

1 

116092 

NRS 

484.2 

222.2 

a 

14  26.47.9 

221.6 

-0  6 

4.1 

3.2 

954  8 

2.56 

1 

116093 

GER 

1077.6 

334.9 

Sn 

14.28  39.2 

2  1 

353  5 

18.6 

4.0 

3.6 

175  5 

1.54 

4 

1 

218 

Origin  time  Lat  Lon  Azres  Timres  Wres  Nphase 


1990 

-292  14 

33  35.0 

48 

40  30.42 

5 

87  0 

22  1  69 

2 

Sta 

Dist 

Az 

Ph 

Time 

Tres 

Azim 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq  Pol 

Arid 

FIN 

1478  5 

167.3 

Pn 

14.36  43  0 

0  0 

155.6 

-11  7 

8.4 

6  5 

94  6 

4  32 

3 

116095 

NRS 

1825  2 

130  7 

Pn 

14  37.24.7 

-0.4 

130  3 

0.0 

10  1 

5.8 

197  2 

4  71 

3 

11609* 

Table  7.3.1.  (Page  2  of  3) 
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U'jawesS'"’ 


j 


Origin  tiae  Lat  Lon  Azres  Timres  Wres  Nphase 

1990-292 -14. 44. 52.0  64.60  37.82  4.45  0.29  1.41  2 


Sta 

Dist 

Az 

Ph 

Tiae 

Tres 

Azin 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq  Pol 

Arid 

ARC 

766.1 

130.0 

Pn 

14.46.32.8 

-0.5 

131.8 

1.8 

8.3 

5.5 

131  2 

8.81 

3 

116109 

ARC 

766.1 

130.0 

Sn 

14  47.52.6 

0.1 

122.9 

-7.1 

4.2 

2.8 

285  2 

8  16 

3  1 

116111 

Origin  tiae  Lat  Lon  Azres  Tiinres  Wres  Nphase 

1990-292:14.56.29.0  49.90  12.11  4.25  0.48  1.54  2 


Sta 

Dist 

Az 

Ph 

Time 

Trea 

Azla 

Ares 

Vel 

Snr 

Amp 

Freq 

Fkq  Pol 

Arid 

GER 

164.9 

316.0 

Pn 

14.56.55.2 

-0.4 

310.5 

-5.5 

6.9 

3.0 

1375.0 

7.64 

3 

221 

GER 

164.9 

316.0 

Lg 

14.57.15.7 

0.6 

319.0 

3.0 

4  2 

7  2 

1416.7 

4  18 

3  -2 

222 

Origin  tiae  Lat  Lon  A zres  T Imres  Wres  Nphase 

1990-292:15. 59.44.0  59.20  11.31  2.60  0.52  1.17  2 


Sta 

Dist 

Az 

Ph 

Time 

Tres 

Azla 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq  Pol 

Arid 

NRS 

171.7 

184.5 

Pn 

16.00.11.8 

0.2 

185.3 

0.8 

7.0 

13  6 

462.6 

4  38 

1 

116126 

NRS 

171.7 

184.5 

So 

16  00.34.5 

-0  8 

180.1 

-4.4 

3.9 

4.7 

864  1 

9.00 

3  -3 

116127 

Origin  tiae  Lat  Lon  Azres  T lares  Wres  Nphase 

1990-292:16.09.58.0  49.60  11.64  5.55  0.29  1.68  2 


Sta 

Dist 

AZ 

Ph 

Tiae 

Tres 

Azla 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq  Pol 

Arid 

GER 

171.9 

300  0 

Pn 

16  10.25  2 

-0.5 

289.5 

-10.5 

6.5 

4.0 

800.7 

9  71 

4 

228 

GER 

171.9 

300.0 

Lg 

16.10.46.0 

-0.1 

300.6 

0.6 

4.0 

3.3 

758  0 

4  95 

3 

229 

Origin  tiae  Lat  Lon  Azres  Timres  Wres  Nphase 

1990-292.20.46  18.0  60  70  21.82  5.40  0.43  1.78  2 


Sta 

Dist 

AZ 

Ph 

Tiae 

Tres 

Aria 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq 

Pol 

Arid 

FIN 

244.7 

252  l 

Pn 

20.46.55  3 

-0.3 

261.7 

9.6 

8  7 

5  7 

115.6 

4  22 

2 

116244 

FIN 

244  7 

252  1 

Sn 

20  47.27.1 

0.5 

253  2 

1.1 

4  6 

3  8 

241  1 

4.66 

3 

-3 

116245 

Origin  tiae  Lat  Lon  Azres  Timres  Wres  Nphase 

1990-292:21.57.40.0  51.70  15.68  4.23  1.70  2.75  6 


Sta 

Dist 

Az 

Ph 

Tiae 

Tres 

Alia 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq 

Pol 

Arid 

GER 

347  5 

23.2 

Pn 

21  58  30.3 

0  0 

29  3 

6.1 

8.9 

64.9 

2858  7 

2  75 

1 

1 

253 

GER 

347.5 

23.2 

P 

21.58.38  2 

25.0 

1  8 

6.4 

17.3 

3017.5 

2  64 

1 

1 

254 

GER 

347  5 

23.2 

P 

21  58  43.8 

25  6 

2  4 

6  3 

8  7 

4004.7 

2  49 

3 

255 

GER 

347  5 

23.2 

Lg 

21  59.25.2 

8  0 

26  2 

3  0 

4.2 

3  4 

8589.2 

2.27 

l 

257 

GER 

347.5 

23.2 

s 

21.59.38.5 

29  8 

6.6 

5.3 

3  2 

4114.4 

1  00 

1 

3 

258 

NRS 

1038.4 

163.9 

Pn 

21.59  54  8 

0  4 

166  3 

2  4 

8  3 

17.1 

383  9 

3.42 

1 

116282 

NRS 

1038.4 

163  9 

P 

22.00.04.8 

162  3 

-1  6 

9.7 

4.5 

245.6 

4.02 

1 

116283 

NRS 

1038.4 

163.9 

Sn 

22  01.38  5 

-0  4 

167  9 

4  0 

4.8 

3  7 

417  7 

2.27 

1 

116284 

FIN 

1256  7 

215.0 

Pn 

22  00  21.3 

0.3 

217.4 

2  4 

9.8 

23.8 

195  6 

3  78 

2 

1 

116281 

FIN 

1256.7 

215  0 

Sn 

22.02  26  0 

1.1 

207  5 

-7  5 

4  9 

2  6 

202  5 

2.84 

1 

116286 

Origin  time  Lat  Lon  Azres  Tiares  Wres  Nphase 


1990- 

292  •  23 

36.38  0 

61 

'.90  19  99 

7. 

23  2 

29  4  09 

5 

Sta 

Dist 

Az 

Ph 

Time 

Tres 

Azla 

Ares 

Vel 

Snr 

Aap 

Freq 

Fkq  Pol 

Arid 

ARC 

288.7 

233.4 

Pn 

23  37.20  6 

-0  4 

225  4 

-8  0 

7.4 

55.2 

1175.6 

6  11 

2  -1 

116311 

ARC 

288.7 

233.4 

P 

23.37.25  5 

221  6 

-11  8 

6.7 

7.3 

695  8 

4.11 

2 

116312 

ARC 

288,7 

233  4 

P 

23  37.46.9 

221.7 

-11  7 

7  1 

2  4 

762  6 

3.67 

1 

116313 

ARC 

288.7 

233  4 

Sn 

23.37.53.4 

-2  5 

227.8 

-5.6 

S  0 

15  2 

2346  3 

3  61 

2 

116314 

ARC 

288.7 

233.4 

23.38.03  4 

4.7 

225  3 

-8.1 

4.0 

3.0 

2834  8 

1  42 

1 

116315 

FIN 

776  3 

340.7 

Lg 

23  40.11.7 

-3.4 

346  2 

5  5 

4.1 

3  0 

198.6 

1  79 

1 

116310 

NRS 

896.3 

23.3 

Pn 

23.38.34.8 

-0.4 

32.3 

9.0 

8.3 

5.4 

114.1 

3  45 

3 

116309 

Origin  time  Lat  Lon  Azres  Tiares  Wres  Nphase 


1990- 

•292  23. 

42  32.0 

67 

.30  20  04 

1. 

54  2.52 

2.90 

3 

Sta 

Dist 

AZ 

Ph 

Tiae 

Tres 

Azla 

A res 

Vel 

Snr 

Aap 

Freq 

Fkq  Pol 

Arid 

ARC 

3:’..  7 

224  6 

Pn 

23.43  21.2 

0  4 

225  2 

0  6 

7.5 

24  7 

979.0 

6  33 

1  -1 

116316 

ARC 

335  7 

224.6 

P 

23.43.23  5 

220  3 

-4  3 

7  3 

6.1 

470  3 

5  04 

2 

116317 

ARC 

335.7 

224  6 

P 

23.43  32  7 

226  8 

2.2 

7  4 

4  8 

443  2 

3  99 

3 

116318 

ARC 

335.7 

224.6 

Sn 

23.43.53.8 

-6.2 

225  5 

0.9 

4  8 

6.3 

944  5 

4  68 

2 

116319 

ARC 

335  7 

224.6 

Lg 

23.44.06.9 

l  0 

221  5 

-3  1 

4.3 

4.8 

1896.3 

2  90 

1 

116320 

Table  7.3.1.  (Page  3  of  3) 
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LATITUDE  (DEG  N) 


Initial  beam  grid 


Fig.  7.3.1.  Initial  beam  grid  used  by  the  generalized  beamforming  method 
when  processing  data  from  four  regional  arrays.  ARCEo  j,  FINESA,  NORESS 
and  GERESS  are  marked  with  solid  triangles. 
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LATITUDE  (DEG  N) 


Events  located  by  two  or  more  arrays 


0.0  10.0  20.0  30.0 

LONGITUDE  (DEG  E) 


Map  projection:  LAMBERT  EQUAL  AREA,  AZIMUTHAL 


Fig.  7.3.2.  Events  located  by  GBF  having  four  or  more  associated  phase 
detections.  The  data  cover  the  five-diy  period  22-26  October  1990.  Out  of 
a  total  of  37  events,  12  were  located  by  two  arrays,  22  by  three  arrays  and  3 
by  four  arrays.  The  four  regional  arrays,  ARCESS,  FINESA,  NORESS  and 
GERESS,  are  marked  with  solid  triangles. 
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7.4  Continuous  threshold  monitoring  of  the  Novaya  Zemlya 
test  site 

Introduction 

The  continuous  threshold  monitoring  technique  (Ringdal  and  Kvaerna, 
1989)  represents  a  new  approach  toward  achieving  reliable  seismic  monitoring 
for  the  purpose  of  verifying  nuclear  test  ban  treaties. 

Traditionally,  seismic  monitoring  has  relied  upon  applying  signal  detec¬ 
tors  to  individual  stations  within  a  monitoring  network,  associating  detected 
phases  and  locating  possible  events  in  the  region  of  interest.  This  procedure 
has  been  accompanied  by  assessments  of  network  capabilities  for  the  target 
region,  usually  by  applying  statistical  models  for  the  noise  level  distribution, 
introducing  station  corrections  for  signal  attenuation  and  devising  a  combi¬ 
national  procedure  to  determine  the  detection  threshold  as  a  function  of  the 
number  of  phase  detections  required  for  reliable  location. 

The  statistical  noise  models  used  in  these  capability  assessments  are  not 
able  to  accommodate  the  effect  of  interfering  signals,  such  as  the  coda  of  large 
earthquakes,  which  may  cause  the  estimated  thresholds  to  be  quite  unrealistic 
at  times.  Furthermore,  only  a  statistical  capability  assessment  is  achieved, 
and  no  indication  is  given  as  to  particular  time  intervals  when  the  possibility 
of  undetected  clandestine  explosions  is  particularly  high. 

The  continuous  threshold  monitoring  technique  alleviates  all  of  these  prob¬ 
lems.  It  makes  it  possible  to  ascertain,  at  any  point  in  time,  for  a  given  target 
region,  the  maximum  magnitude  of  a  possible  clandestine  explosion  at  a  pre¬ 
defined  level  of  confidence.  This  makes  it  oossible  to  focus  attention  upon 
those  specific  time  intervals  when  realistic  t  'asion  opportunities  exist,  while 
retaining  confidence  that  no  treaty  violation  has  occurred  at  other  times. 

Application  to  the  Novaya  Zemlya  test  site 

In  order  to  demonstrate  how  such  monitoring  could  be  performed  in  a 
practical  operational  situation,  we  have  conducted  an  experiment  during  which 
we  have  applied  continuous  threshold  monitoring  to  the  Novaya  Zemlya  test 
site  for  a  full  one-week  period.  Our  data  base  has  been  the  Fennoscandian 
regional  array  network  (NORESS,  ARCESS,  FINESA).  As  illustrated  in  Fig. 
7.4.1,  these  three  arrays  are  all  within  regional  distances  from  the  test  site,  with 
excellent  P-phase  detection  capabilities  (Fig.  7.4.2).  The  ARCESS  array  also 
detects  S  phases  from  Novaya  Zemlya  explosions  quite  well,  whereas  NORESS 
and  FINESA  have  a  lower  S-phase  detection  capability. 

The  parameters  used  in  the  threshold  monitoring  experiment  are  given  in 
Table  7.4.1.  For  each  array,  we  steer  “optimum”  P  and  S  beams  towards  the 
test  site,  and  calibrate  these  beams  using  actually  observed  signal  attenuation 
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from  Novaya  Zemlya  explosions.  By  focusing  in  this  way  on  the  target  region, 
we  can  at  any  point  in  time  measure  the  “noise  magnitude”  for  a  given  phase 
at  a  given  array,  and  combine  these  data  to  obtain  a  network  threshold  as 
explained  in  detail  by  Ringdal  and  Kvaerna  (1989). 

Results 

Figs.  7.4. 3-7.4. 9  show  the  result  of  the  monitoring  experiment.  Each  of 
these  figures  covers  one  data  day,  starting  24  October  1990.  The  upper  three 
traces  of  each  figure  represent  the  thresholds  (i.e.,  90  %  upper  magnitude 
limits)  obtained  from  the  three  individual  arrays,  whereas  the  bottom  trace 
illustrates  the  network  threshold.  Typically,  the  individual  array  traces  have  a 
number  of  significant  peaks  for  each  24-hour  period,  due  to  interfering  events 
(local  or  teleseismic).  On  the  network  trace,  the  number  and  sizes  of  these 
peaks  are  greatly  reduced,  because  an  interfering  event  will  usually  not  provide 
matching  signals  at  all  the  stations.  From  probabilistic  considerations,  it  can 
in  such  cases  be  inferred  that  the  actual  network  threshold  is  lower  than  these 
individual  peaks  might  indicate. 

On  each  of  the  one-day  figures,  we  have  included  comments  explaining  the 
presence  of  the  most  significant  peaks  on  the  network  trace.  Here,  we  will 
just  note  that  the  first  day,  24  October  1990,  was  the  day  of  an  actual  nuclear 
explosion  (mt  =  5.7)  on  Novaya  Zemlya,  and  this  event  naturally  stands  out 
on  the  plot.  While  the  peak  value  of  the  network  threshold  plot  does  not 
represent  the  actual  magnitude  of  the  event,  it  is  in  fact  quite  dose  (5.64). 

As  a  general  comment  to  Figs.  7.4. 3-7.4. 9,  we  note  that  such  plots,  which 
are  easily  generated  by  the  Intelligent  Monitoring  System  (IMS)  (Bache  et  al, 
1990),  will  enable  the  analyst  to  obtain  an  instant  assessment  of  the  actual 
threshold  level  of  the  monitoring  network.  The  peaks  on  the  network  traces 
may  be  quickly  correlated  with  the  IMS  detection  bulletin,  in  order  to  decide 
whether  they  originate  from  interfering  events  or  from  events  in  the  target 
region. 

Discussion 

In  a  monitoring  situation,  it  will  be  important  to  isolate  and  analyze  more 
extensively  those  time  intervals  which  offer  significant  evasion  opportunities. 
Table  7.4.2  gives  a  statistic  of  the  number  of  occasions  during  which  the  up¬ 
per  magnitude  limit  exceeded  a  given  level.  In  theory,  if  this  limit  is,  e.g.,  at 
3.0,  it  might  be  possible  that  a  clandestine  m &  =  3.0  explosion  had  occurred 
without  being  detected.  There  are  many  options  available  to  investigate  such 
a  hypothesis  in  more  detail,  although  we  have  not  attempted  to  do  so  in  this 
study.  The  most  immediate  approach  would  be  to  analyze  high-frequency  sig¬ 
nals  for  the  time  interval  being  considered.  For  example,  cn  ARCESS  records 
Novaya  Zemlya  explosions  will  contain  significant  energy  at  10  Hz  and  above, 
even  at  magnitudes  well  below  3.0.  Teleseismic  events,  even  of  large  mj,,  will 


83 


not  contain  much  energy  at  these  frequencies  and  thus  it  might  be  possible  to 
obtain  additional  indications  from  these  data. 

To  assess  interfering  phases  from  events  at  regional  distances  is  more  diffi¬ 
cult,  since  the  high-frequency  energy  might  not  discriminate  such  events  from 
Novaya  Zemlya  explosions.  In  such  cases,  additional  procedures,  such  as  max¬ 
imum  likelihood  beamforming,  might  become  useful  to  suppress  signals  from 
the  interfering  event  and  thereby  obtain  a  more  realistic  estimate  of  the  signal 
energy  arriving  from  the  target  region. 

It  is  significant  that  the  3-array  network  studied  in  this  paper  can  monitor 
the  Novaya  Zemlya  test  site  down  to  m b  2.5  or  below  more  than  99  %  of  the 
time  (Fig.  7.4.10).  Further  improvements  would  clearly  be  possible  by  adding 
more  stations  to  the  monitoring  network,  especially  highly  sensitive  stations 
at  other  azimuths  than  those  covered  by  the  Fennoscandian  network.  This 
would  in  particular  contittBite  to  lowering  the  peaks  due  to  interfering  events, 
whereas  any  event  truly  originating  in  the  target  region  would  of  course  still 
stand  out  clearly  on  the  combined  network  traces. 

In  conclusion,  the  continuous  threshold  monitoring  has  been  demonstrated 
to  provide  a  simple  and  very  effective  tool  in  day-to-day  monitoring  of  a  site 
of  particular  interest.  Further  research  will  focus  upon  developing  methods 
to  analyze  time  intervals  during  which  significant  evasion  possiblities  might 
exist.  Data  from  the  regional  arrays,  the  large-aperture  NORSAR  arrray,  as 
well  as  other  available  stations,  will  be  used  in  these  analyses. 


T.  Kvaerna 
F.  Ringdal 


References 

Bache,  T.C.,  S.R.  Bratt,  J.  Wang,  R.M.  Fung,  C.  Kobryn  and  J.  Given  (1990): 
The  Intelligent  Monitoring  System,  Bull.  Seism.  Soc.  Am.,  Special  Issue, 
in  print. 

Ringdal,  F.  and  T.  Kvaerna  (1989):  A  multichannel  processing  approach  to 
real  time  network  detection,  phase  association  and  threshold  monitoring, 
Bull.  Seism.  Soc.  Am.,  79,  1927-1940. 


84 


V1 


* 


Station 

Phase 

Tr.  Time 

App.  Vel. 

Azim. 

Filter 

Config. 

STAJen. 

Tim.  Tol. 

Sta.Calib 

ARC 

Pn 

148.0 

9.9 

60.5 

3.0-5.0 

A0,B,C,D 

2.0 

2.0 

0.754 

ARC 

Sn 

257.0 

4.9 

53.2 

3.0-5. 0 

A0,B,C,D 

5.0 

3.0 

1.176 

FIN 

P 

228.0 

9.6 

32.9 

2.0-4.0 

A0,B,C, 

2.0 

2.0 

1.520 

NRS 

P 

284.0 

10.4 

28.1 

1.5-3. 5 

A0,B)C,D 

2.0 

2.0 

0.677 

Tr.  time 
App.  vel. 
Azim. 
Filter 
Config. 

STAJen. 
Tim.  tol. 
STA_calib 


Travel  time  of  phase 

Apparent  velocity  from  broadband  F-K  measurement 
Azimuth  from  broadband  F-K  measurement 
Cutoffs  of  bandpass  filter  (3rd  order  Butterworth) 

Array  configuration  used  in  beamforming.  A0,B,C  means 
AOZ,  B-ring  and  C-ring 
STA  length  in  seconds 

Time  tolerance  when  searching  for  maximum  STA 
Calibration  factor  used  when  converting  STA  values 
(in  quantum  units)  to  magnitude. 

Magnitude  =  loglO(STA)  -f  STA.calib. 


Table  7.4.1. 


297 

298 

Day-of-Year 
299  300  301 

302 

303 

Total 

(one  week) 

mfc  >  4.0 

1 

0 

0 

0 

0 

0 

0 

1 

mt  >  3.5 

1 

1 

0 

0 

0 

0 

0 

2 

mi  >  3.0 

3 

2 

2 

0 

0 

0 

0 

7 

mi  >  2.5 

5 

12 

5 

0 

3 

6 

3 

34 

Table  7.4.2.  Statistics  of  peaks  in  the  network  threshold  traces.  See  also  '•omments  on 
Figs.  7.4.3-7.4.9. 
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LATITUDE  (DEG  N) 


Novaya  Zemlya 


Fig.  7.4.1.  Location  of  the  target  area  (Novaya  Zemlya)  for  the  threshold 
monitoring  experiment.  The  locations  of  the  three  arrays  NORESS  (A  =  2280 
km),  ARCESS  (A  =  1110  km)  and  FINESA  (A  =  1780  km)  are  indicated. 


86 


Noise  x  1000 


.  1  5E-+-05I 


ARC 


1A.99.S2.-lBO 


3-5  HZ 


37S97 


FIN 


2.05E1-0SP 


NR 


Noise  x  100 


Noise  x  100 

vv4aa'V\'vv/’W'~YA/vvvw.^w-'V^ — 


2-4  Hz 


19.02.07.010 


1. 5-3.5  Hz 


0.000 


20.000 


40.000 


60.000  Sec 


1 990-297 


5/90  11:10:38  NORSAR 


Fig.  7.4.2.  P  -wave  recordings  (filtered  array  beams)  at  ARCESS,  FINESA 
and  NORESS  for  the  Novaya  Zemlya  nuclear  explosion  24  October  1990.  To 
illustrate  the  high  signal-to-noise  ratios,  noise  traces  scaled  by  factors  of  1000 
(ARCESS)  and  100  (FINESA  and  NORESS)  are  displayed  for  each  array. 
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Fig.  7.4.3.  Threshold  monitoring  of  the  Novaya  Zemlya  test  site  for  day 
297  (24  October  1990).  The  top  three  traces  represent  thresholds  (upper  90 
per  cent  magnitude  limits)  obtained  from  each  of  the  three  arrays  (AECESS, 

FINESA,  NORESS),  whereas  the  bottom  trace  shows  the  combined  network 
thresholds. 

Notes: 

1.  An  underground  nuclear  explosion  (mb  =  5.7)  at  Novaya  Zemlya 
at  14.58.00  GMT.  The  peak  of  the  network  trace  is  5.64. 

2.  Two  teleseismic  earthquakes  from  N.  Xinjang  province,  China 
(mj  =  5.2  and  5.4).  The  P-wave  train  from  each  of  these 
earthquakes  causes  the  network  threshold  to  increase  to  about 
m*  =  3.0  for  the  target  region.. 
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Fig.  7.4.4.  Same  as  Fig.  7.4.3,  but  for  day  298  (25  October  19p0).  The 
FINESA  array  had  several  short  outages  this  day,  but  this  caused  no  particular 
problems  in  terms  of  network  threshold  capacity. 

Notes: 

3.  An  earthquake  (m;,  =  4.5)  near  Jan  Mayen.  The  corresponding 
network  threshold  peak  fo.-  Novaya  Zemlya  is  mf,=2.8. 

4.  A  teleseismic  earthquake  (m;,  -  6.0)  at  Hindu  Kush.  The  relatively 
strong  P-wave  train  caused  a  peak  threshold  of  mj  =  3.8  for  monitoring 
Novaya  Zemlya. 

5.  A  teleseismic  earthquake  (mj  =  5.9)  at  Mindanao,  Philippine 
Islands.  Corresponding  threshold  is  mj  =  3.0. 

6-7.  A  sequence  of  seismic  events  (presumably  underwater  explosions) 

near  Murmansk,  Kola  Peninsula.  The  network  threshold  for  monitoring 
Novaya  Zemlya  is  about  mt  =  2.5  to  2.8  at  the  times  of  these 
events. 
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Fig.  7.4.5.  Same  as  Fig.  7.4.3,  but  for  day  299  (26  October  1990). 

Notes: 

8.  A  large  mining  explosion  (Ml  =  3.0)  near  the  Norway-USSR  border. 

The  Novaya  Zemlya  threshold  peak  is  m;,  —  3.2.  Note  that  this  threshold 
exceeds  the  event  magnitude;  this  is  beca  use  of  the  proximity  of  the  event 
to  the  network  stations. 

9.  A  large  mining  explosion  (Ml  =  2.5)  at  the  Kola  Peninsula.  The 
Novaya  Zemlya  threshold  peak  is  m&  -  2.8. 

10.  A  teleseismic  earthquake  (mj  =  5.1)  near  Lake  Baikal.  The  network 
threshold  peak  is  mj  =  3.2. 
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Fig.  7.4.6.  Same  as  Fig.  7.4.3,  but  for  day  300  (27  October  1900). 

Notes: 

No  significant  peak  during  this  day.  The  simultaneous  outages  at  FINESA 
and  NORESS  did  not  cause  significant  problems  for  the  network  threshold 
monitoring. 


Magnitude 


Fig.  7.4.10.  Cumulative  statistics  of  the  network  threshold  magnitudes  from 
Figs.  7.4.3-7.4.9,  covering  one  full  week  of  data. 


7.5  Real-time  processing  using  a  hybrid  3-component  /  small 
array  station 

Introduction 

In  recent  years,  much  effort  has  been  devoted  to  developing  and  assessing 
various  methods  for  3-component  seismic  data  processing.  It  has  been  demon¬ 
strated  that  polarization  analysis  can  provide  P  azimuth  estimates  v/ith  good 
accuracy  from  a  single  3-component  station.  Using  SH  or  SV  particle  motion 
models,  some  success  has  also  been  reported  in  determining  azimuth  from  S 
and  Lg  phases,  although  there  is  often  a  90  or  180  degree  ambiguity  in  the 
resulting  estimates. 

These  efforts  notwithstanding,  a  fundamental  problem  in  automatic  3- 
component  processing  still  remains.  This  is  the  problem  of  phase  identification, 
which  of  course  is  essential  in  deciding  which  particle  motion  model  to  apply. 
Prom  our  experience  at  NORESS,  a  high  degree  of  rectilinearity,  which  in 
theory  would  indicate  the  presence  of  a  P  phase,  is  quite  often  seen  also  for  S 
and  Lg  phases,  and  even  for  noise  bursts,  when  applying  polarization  analysis. 
As  a  result,  an  automatic  P-type  solution  found  by  3-component  processing 
is  often  of  limited  usefulness,  since  it  may  either  be  a  P  phase  (with  correct 
azimuth)  or  some  other  phase  (in  which  case  the  azimuth  estimate  is  useless). 

Small-aperture  arrays  of  the  NORESS-type  have  proved  to  be  very  effective 
in  processing  regional  as  well  as  teleseismic  phases.  Their  primary  features  are: 

-  Significant  SNR  gains  at  high  frequencies 

-  Reliable  phase  identification  (P-type  versus  S-type  phases) 

-  Precise  azimuth  estimates  of  all  phase  types. 

While  the  accuracy  of  NORESS  azimuth  estimates  can  be  as  good  as  ± 
1  degrees  for  well-calibrated  regions,  the  uncertainty  of  uncalibrated  azimuth 
estimates  is  often  of  the  order  of  10  degrees  or  more,  due  to  lateral  inho¬ 
mogeneities  near  the  receivers.  In  practical  schemes  for  automatic  phase  as¬ 
sociation  (e.g.,  Mykkeltveit  and  Bungum,  1984),  a  tolerance  of  30  degrees  in 
azimuth  deviation  from  the  true  value  is  often  assumed.  Given  that  the  accept¬ 
able  limits  in  azimuth  estimation  for  phase  association  purposes  are  much  less 
restrictive  than  the  optimum  array  capability,  a  natural  question  is  whether 
a  smaller  regional  array  can  achieve  reliable  phase  identification  as  well  as  an 
acceptable  uncertainty  in  azimuth  estimation. 

We  have  attempted  to  address  this  question,  and  have  chosen  to  evaluate 
the  smallest  such  array  available  to  us:  the  NORESS  A-ring  subarray  (Fig. 
7.5.1).  This  subarray  comprises  a  center  3-component  seismometer  AO,  sur¬ 
rounded  by  a  triangular  SPZ  array  Al-3.  The  diameter  of  this  subarray  is 
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only  0.3  kr  ,  i.e.,  a  factor  of  10  less  than  NORESS,  and  it  thus  spans  an  area 
only  1  per  cent  of  that  of  the  full  NORESS  array.  In  this  paper,  we  present 
an  initial  evaluation  of  automatic  detection  processing  using  this  very  small 
array. 

Automatic  processing 

We  conducted  automatic  detection  processing  of  the  A-ring  subarray  for 
a  period  of  5  days  (22  -26  October  1990),  using  a  standard  multi-filter  power 
detector  with  parameters  as  specified  in  Table  7.5.1.  For  each  detected  signal, 
broadband  F-K  analysis  was  conducted  using  the  4  SPZ  sensors,  and  compared 
to  the  NORESS  F-K  results.  We  also  carried  out  polarization  analysis  (using  a 
P-wave  particle  motion  model)  for  each  detected  phase,  applying  the  method 
of  Kvaerna  and  Doornbos  (1986)  to  the  3-component  AO  system. 

To  obtain  a  data  base  for  the  evaluation,  we  extracted  all  seismic  phases 
detected  by  NORESS  and  associated  to  verified  regional  events  for  the  5-day 
period.  The  generalized  beamforming  procedure  (Ringdal  and  Kvaerna,  1989) 
and  the  results  from  IMS  processing  (Bache  et  al,  1990)  were  used  in  order 
to  ensure  correct  location  of  these  reference  events.  P-coda  detections  and 
multiple  S-phases  were  ignored,  so  that  each  event  provided  a  maximum  of  3 
phases  (P,  S  and  Lg).  These  phases  were  then  matched  to  the  detection  lists 
produced  from  the  A-ring  SPZ  subarray  and  the  3-component  analysis,  and 
the  phase  velocity  and  azimuth  estimates  were  compared. 

Phase  identification 

Fig.  7.5.2  shows  the  A-ring  phase  velocity  F-K  estimates  for  P  phases 
(circles)  and  S  phases  (crosses)  for  the  reference  data  set.  The  separation  is 
better  than  95  per  cent,  which  implies  that  even  this  small  array  is  able  to 
provide  correct  phase  identification  automatically  and  with  high  confidence. 
The  corresponding  picture  for  the  full  NORESS  array  is  shown  in  Fig.  7.5.3, 
and  naturally  shows  a  better  separation.  It  should  be  noted  here  that  all  the 
reference  phases  were  by  definition  correctly  identified  by  NORESS,  so  it  is 
not  surprising  that  the  separation  is  100  per  cent.  But  it  is  still  interesting  to 
compare  the  scatter  in  the  two  data  sets. 

P-wave  azimuths 

Fig.  7.5.4  compares  P-wave  azimuths  estimated  by  the  full  NORESS  SPZ 
array  and  the  small  A-ring  SPZ  subarray  using  broadband  F-K  analysis  in 
both  cases.  The  estimates  are  very  consistent,  and  the  large  majority  are  well 
within  a  tolerance  limit  of  30  degrees. 

A  corresponding  plot  for  P  waves  analyzed  from  3-component  records  is 
given  in  Fig.  7.5.5,  and  shows  a  similar  amount  of  scatter.  In  fact,  the  scatter 
is  slightly  less  than  in  Fig.  7.5.4.  Thus,  the  polarization  analysis  can  give 
very  useful  contributions  to  regional  azimuth  estimates,  assuming  that  the 
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phase  first  has  been  identified  as  a  P  phase.  An  interesting  approach  would 
be  to  combine  the  two  azimuth  estimates  (A-ring  SPZ  and  3-comp)  in  order 
to  obtain  increased  stability. 

S-wave  azimuths 

Fig.  7.5.6  compares  azimuths  of  S-type  phases  estimated  by  the  A-ring  SPZ 
subarray  and  the  full  NORESS  array.  Again,  the  correspondence  is  quite  good, 
and  thus  shows  that  automatic  regional  phase  association  using  RONAPP- 
type  processing  can  be  achieved  even  when  using  an  array  as  small  as  0.3 
km  in  diameter.  We  note  that  our  3-component  processing  provides  no  useful 
azimuth  information  for  S  phases,  but  it  is  of  course  possible  that  such  infor¬ 
mation  could  be  obtained  in  certain  cases,  given  that  the  phase  first  has  been 
identified  as  S  or  Lg. 

Detectability 

Figs.  7.5.7  and  7.5.8  illustrate  the  P  and  S  wave  detectability  of  the  A- 
ring  array  as  a  function  of  NORESS  SNR.  From  Fig.  7.5.7,  it  is  seen  that 
all  P  phases  with  SNR  >  20  dB  (i.e.,  STA/LTA  >  10  at  NORESS)  have 
been  detected.  At  distances  below  500  km,  several  events  of  relatively  low 
SNR  at  NORESS  have  also  A-ring  detections;  this  is  due  to  the  high  signal 
frequencies  which  cause  the  full  array  SNR  gains  of  these  phases  to  be  less 
than  the  theoretical  y/N.  At  distances  above  500  km,  the  superiority  of  the 
full  NORESS  array  becomes  apparent. 

In  Fig.  7.5.8  it  is  seen  that  the  A-ring  subarray  is  close  to  matching 
NORESS  S-phase  detectability  at  all  distances.  This  is  because  the  horizontal 
components  of  the  A0  3-component  system  provide  quite  efficient  detection  of 
S  and  Lg  phases,  in  particular  when  added  incoherently  to  the  vertical  com¬ 
ponent.  The  full  array  does  not  have  the  same  SNR  gain  for  secondary  phases 
as  for  P-phases,  because  of  less  signal  coherency  and  (in  particular)  coherent 
noise  caused  by  P  coda.  Thus,  as  expected,  there  is  not  a  significant  difference 
in  secondary  phase  detection. 

False  alarm  consideration 

In  practical  operation  of  any  seismic  system,  the  problem  of  false  detections 
is  very  important.  This  is  especially  the  case  if  the  online  detector  is  operated 
at  a  low  detection  threshold,  and  it  is  essential  to  be  able  to  identify  false 
alarms  at  as  early  a  stage  as  possible. 

To  address  this  problem,  we  have  analyzed  in  detail  all  the  A-ring  detec¬ 
tions  for  one  full  data  day  (24  October  1990).  The  results  are  presented  in 
Table  7.5.2,  again  with  NORESS  results  as  a  reference.  From  the  table,  it  is 
seen  that  out  of  117  total  phase  detections,  105,  or  90  per  cent,  were  correctly 
classified  using  the  broadband  F-K  analysis  applied  to  A-ring  SPZ  data.  Of 
these  105  phases,  37  were  P,  36  S  (or  Lg)  and  32  noise  (i.e.,  low-velocity  de- 
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tections).  Note  that  P  coda  detections  were  counted  as  P  and  S  v.oda  were 
counted  as  S  in  these  statistics.  None  of  the  37  phases  which  (according  to 
NORESS)  were  of  the  P  type  were  misclassified  by  the  A-ring.  Of  the  40  S 
phases,  4  were  misclassified  (2  P  and  2  noise).  Out  of  40  noise  detections,  3 
were  given  P-phase  velocities  and  5  were  given  S-phase  velocities  when  using 
the  A-ring. 

These  statistics  must  be  considered  quite  satisfactory.  In  fact,  it  appears 
that  the  SNR  threshold  for  the  A-ring  detector  could  be  lowered  significantly, 
and  still  produce  a  reasonable  false  alarm  rate. 

Conclusions 

The  problems  encountered  when  using  a  3-component  system  in  a  real¬ 
time  automatic  processing  environment  appear  to  be  effectively  alleviated  by 
supplementing  the  3-C  system  with  a  very  small  (aperture  0.3  km)  3-element 
array.  Based  on  the  studies  in  this  paper  of  the  NORESS  A-ring  subarray,  we 
conclude  that: 

-  Phase  identification  (P  or  S)  can  be  reliably  achieved  (better  than  95  % 
success  rate)  using  F-K  analysis  of  the  A-ring  array  data. 

-  Azimuth  estimates,  with  accuracy  generally  to  within  30  degrees  can  be 
obtained  for  both  P  and  S  phases  using  A-ring  array  processing,  and  (at 
least)  for  P-phases  using  3-component  polarization  analysis  of  the  AO 
3-component  station.. 

-  Good  regional  P-phase  detectability  can  be  obtained  from  the  A-ring 
array  out  to  500  km  epicentral  distance.  At  greater  distances,  P-wave 
detectability  relative  to  that  of  NORESS  deteriorates  sharply. 

-  Detectability  of  S-phases  using  the  A-ring  array  is  excellent  at  all  dis¬ 
tances,  and  comes  close  to  matching  that  of  the  full  NORESS  array. 

-  The  A-ring  array  analysis  makes  it  possible  to  isolate  the  majority  of 
noise  detections,  thus  giving  an  acceptable  false  alarm  rate  for  online 
operation. 

It  would  be  an  interesting  experiment  to  supplement  some  existing  3- 
component  stations  with  a  triangular  small  array  as  described  here,  and  con¬ 
duct  network  detection  and  location  experiments.  Before  doing  so,  it  is  of 
course  important  to  verify  that  these  results  can  be  obtained  in  other  geolog¬ 
ical  and  geographical  environments,  e.g.,  by  analyzing  similar  data  for  other 
existing  arrays  (ARCESS,  GERESS,  FINES  A). 
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A  very  smaii  euray  of  the  A-ring  type  is  especially  suited  for  processing 
high  signal  frequencies.  In  fact,  one  might  consider  a  dense  sensor  deploy¬ 
ment  within  the  A-ring  aperture,  as  an  experiment  aimed  at  conducting  array 
processing  at  frequencies  of  20-40  Hz.  This  would  of  course  require  a  higher 
sampling  rate  than  ihe  40  Hz  currently  used  at  NORESS.  Array  processing 
at  these  frequencies  would  be  of  particular  interest  in  the  context  of  devel¬ 
oping  methods  for  monitoring  cavity  decoupled  explosions,  which  might  have 
significant  signal  energy  in  this  frequency  band. 


T.  Kvaerna 
F.  Ringdal 
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Coherent  beams: 


Beam 

Apparent  vel. 

Azimut 

Filter 

band 

Configuration 

Threshold 

NZ01 

9999.9 

0.0 

0.5- 

1.5 

AOz, Alz, A2z, A3z 

4.0 

NZ02 

9999.9 

0.0 

1.0- 

2.0 

A0z,Alz,A2z,A3z 

3.7 

NZ03 

9999.9 

0.0 

1.0- 

3.0 

AGz,Alz,A2z,A3z 

3.7 

NZ04 

9999.9 

0.0 

1.5- 

2.5 

AOz , Alz , A2z , A3z 

3.5 

NZ05 

9999.9 

0.0 

1.5- 

3.5 

A0z,Alz,A2z,A3z 

3.5 

NZ06 

9999.9 

0.0 

2.0- 

3.0 

A0z,Alz,A2z,A3z 

3.5 

NZ07 

9999.9 

0.0 

2.0- 

4.0 

AOz , Alz , A2z , A3z 

3.5 

NZ08 

9999.9 

0.0 

2.5- 

4.5 

AOz, Alz, A2z, A3z 

3.5 

NZ09 

9999.9 

0.0 

3.0- 

5.0 

A0z,Alz,A2z,A3z 

3.5 

NZ10 

9999.9 

0.0 

3.5- 

5.5 

A0z,Alz,A2z,A3z 

3.7 

NZ11 

9999.9 

0.0 

4.0- 

8.0 

AOz ,  Alz ,  A2z ,  A.3z 

3.7 

NZ12 

9999.9 

0.0 

5.0-: 

10.0 

A0z,Alz,A2z,A3z 

3.7 

NZ13 

9999.9 

0.0 

8.0-: 

16.0 

AOz , Alz , A2z , A3z 

4.0 

Incoherent  beams: 


Beam 

Apparent  vel. 

Azimut 

Filter  band 

Configuration 

Threshold 

NA01 

9999.9 

0.0 

0.5-  1.5 

A0z,A0n, AOe 

2.7 

NA02 

9999.9 

0.0 

1.0-  2.0 

A0z,A0n, AOe 

2.6 

NA03 

9999.9 

0.0 

1.0-  3.0 

A0z,A0n,A0e 

2.5 

NA04 

9999.9 

0.0 

1.5-  2.5 

A0z,A0n,A0e 

2.5 

NA05 

9999.9 

0.0 

1.5-  3.5 

A0z,A0n, AOe 

2.5 

NA06 

9999.9 

0.0 

2.0-  3.0 

A0z,A0n,A0e 

2.6 

NA07 

9999.9 

0.0 

2.0-  4.0 

A0z,A0n,A0e 

2.8 

NA08 

9999.9 

0.0 

2.5-  4.5 

A0z,A0n, AOe 

3.4 

NA09 

9999.9 

0.0 

3.0-  5.0 

A0z,A0n,A0e 

3.5 

NA10 

9999.9 

0.0 

3.5-  5.5 

A0z,A0n,A0e 

2.8 

NA1 1 

9999.9 

0.0 

4.0-  8.0 

A0z,A0n,A0e 

2.5 

NA12 

9999.9 

0.0 

5.0-10.0 

A0z,A0n, AOe 

2.5 

NA13 

9999.9 

0.0 

8.0-16.0 

A0z,A0n, AOe 

2.8 

Table  7.5.1.  Parameters  used  for  the  A-ring  detector  experiment. 


Classified  as: 

Correct 

P 

S  or  Lg 

Noise 

Phase  id 

(v>6  km/s) 

(3.4<v<6  km/s) 

(v<3.4  km/s) 

P 

37 

0 

0 

S  or  Lg 

2 

36 

2 

Noise 

3 

5 

32 

Total  phases  detected  by  the  A-ring  detector  :  117 
Total  phases  correctly  classified  :  105  (90  %) 


Table  7.5.2.  Statistics  of  detected  phases  for  the  A-ring  array  for  a  24-hour 
period.  The  phases  are  classified  based  on  estimated  phase  velocities  using  A- 
ring  SPZ  broadband  F-K,  and  the  “correct”  phase  id  is  based  on  F-K  results 
from  the  full  NORESS  array. 
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Fig.  7.5.1.  Geometry  of  the  full  NORESS 
used  in  this  study. 
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Fig,  7,5.2.  Estimated  phase  velocities  using  the  A-ring  SPZ  array  (broadband 
F-K)  for  detected  P  phases  (circles)  and  S  phases  (asterisks).  Note  that  the 
phases  can  be  identified  from  phase  velocity  with  more  than  95  %  accuracy. 
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Fig.  7.5.3.  Similar  to  Fig.  7.5.2,  but  with  phase  velocities  estimated  using 
the  full  NORESS  array. 
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Fig.  7.5.4.  Comparison  of  estimated  azimuths  of  P  phases  using  the  full 
NORESS  array  and  the  A-ring  SPZ  array  (broadband  F-K).  Note  the  excellent 
consistency. 
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Fig.  7.5.5.  Comparison  of  estimated  azimuths  of  P  phases  using  the  full 
NORESS  array  (broadband  F-K)  and  the  AO  3-component  system  (polariza¬ 
tion  analysis).  Note  that  the  consistency  is  similar  to  that  of  Fig.  7.5.4. 
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Fig.  7.5.6.  NORESS  and  A-ring  azimuth  comparison  for  S  phases.  Note  that 
the  consistency  is  as  good  as  for  P  phases  (in  Fig.  7.5.4). 


Fig.  7.5.7.  Illustration  of  P-phase  detectability  of  the  A-ring  SPZ  array 
combined  with  the  3-component  AO  system.  P-phases  detected  by  the  A- 
ring  SPZ/3-component  system  are  marked  as  asterisks,  whereas  nondetected 
phases  are  marked  as  circles.  Note  that  the  reference  array  (NORESS)  is 
clearly  superior  at  distances  >  500  km. 
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|  Fig.  7*5,8*  Same  as  Fig.  7.5.7,  but  for  S  phases.  Note  that  in  this  case 

|  the  sma11  array /3- component  system  comes  close  to  matching  the  full  array 
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7.6  Development  of  two  three-component  stations  in  Poland 
Background 

A  network  of  four  regional  arrays  currently  contributes  data  to  the  NOR- 
SAR  Data  Processing  Center  in  Norway.  These  arrays  are  the  NORESS  and 
ARCESS  arrays  in  Norway,  the  FINESA  array  in  Finland  and  the  GERESS 
array  in  the  Federal  Republic  of  Germany  (see  Fig.  7.6.1).  Data  from  ARCESS 
and  NORESS  are  processed  jointly  in  the  Intelligent  Monitoring  System  (IMS) 
that  has  been  installed  at  NORSAR  this  year.  A  forthcoming  upgrade  of  IMS 
will  allow  joint  processing  of  data  from  N  arrays  and  M  single  stations. 

As  a  cooperative  effort  between  Poland  and  the  United  States,  two  mod¬ 
ern  seismic  stations,  including  workstations  for  on-site  data  recording  and 
analysis,  are  now  being  installed  in  Poland.  Each  installation  comprises  a 
three-component  broad  band  station  and  a  three-component  short  period  sta¬ 
tion.  These  stations  represent  a  valuable  extension  of  the  existing  network  of 
regional  arrays  in  northeastern  Europe.  In  order  for  IMS,  however,  to  take  full 
advantage  of  these  new  stations,  the  data  from  Poland  must  be  transmitted 
continuously  and  in  real  time  to  Norway.  This  contribution  gives  descriptions 
of  the  field  systems  installed  in  Poland,  the  communications  arrangements  be¬ 
tween  Poland  and  Norway,  and  the  system  configured  at  NORSAR  for  the 
purpose  of  acquiring  the  data  from  the  two  stations  in  Poland. 

Installations  in  Poland 

The  two  stations  in  Poland  are  located  at  Ksiaz,  in  the  vicinity  of  the  town 
of  Walbrzych  in  southwestern  Poland,  and  at  Stary  Folwark,  near  the  town  of 
Augustow  in  northeastern  Poland  (see  Fig.  7.6.1). 

Teledyne  Geotech  installed  in  June  this  year  a  seismic  system  at  Ksiaz, 
comprising  a  stiort  period  three- component  station  (GS-13  sensors),  a  broad 
band  three-component  station  (BB-13  sensors),  an  RDAS-200  data  acquisition 
unit,  and  a  SUN-2  based  NOMAD  workstation  for  local  data  analysis  and 
archiving.  Ksiaz  is  an  existing  seismologica)  observatory  affiliated  with  the 
Department  of  Seismology  of  the  Institute  of  Geophysics  (in  Warsaw)  of  the 
Polish  Academy  of  Sciences.  The  local  staff  at  the  observatory  at  Ksiaz  will 
operate  the  new  station  and  use  the  workstation  in  their  daily  data  analysis. 

The  new  station  in  Stary  Folwark  will  be  installed  by  Teledyne  in  Novem¬ 
ber.  It  will  be  identical  to  the  station  at  Ksiaz,  except  there  will  be  no  work¬ 
station,  since  there  will  be  no  local  seismological  staff. 

Communications  links 

Reliable  communications  links  between  the  NORSAR  Data  Processing 
Center  in  Norway  and  the  two  sites  in  Poland  can  only  be  accomplished  via 
satellite.  Discussions  were  conducted  with  the  telecommunication  authorities 


in  both  Poland  and  Noway  in  order  to  determine  how  such  links  could  be  es¬ 
tablished,  and  the  Polish  PTT  responded  that  they  would  allow  the  Norwegian 
Telecommunication  Administration  (NTA)  to  deliver  and  put  into  operation 
the  ground  station  equipment  needed  in  Poland. 

NTA  decided  to  use  the  EUTELSAT  SMS  system  for  the  satellite  links 
between  Poland  and  Norway,  and  communications  were  installed  in  Octo¬ 
ber/November  as  depicted  in  Fig.  7.6.2.  The  data  transmission  links  each 
have  a  capacity  of  64  kbits/s.  The  link  between  Stary  Folwark  and  NORSAR 
is  simplex,  while  the  Ksiaz-NORSAR  link  is  full  duplex.  The  backlink  from 
NORSAR  to  Ksiav:  provides  the  staff  at  the  Ksiaz  observatory  with  access  to 
data  bases  at  NORSAR,  including  Stary  Folwark  data.  The  Ksiaz  to  NOR¬ 
SAR  link  uses  a  set  of  multiplexers  to  enable  division  of  the  capacity  into 
several  channels,  as  shown  in  Fig.  7.6.2.  The  seismic  data  from  Ksiaz  use  a 
19.2  kbits/s  channel,  the  other  19.2  kbits/s  channel  will  be  running  the  SLIP 
protocol  and  used  for  general  purposes,  one  of  the  9.6  kbits/s  channels  will 
be  used  lor  general  terminal  access,  and  two  remaining  channels  (9.6  and  2.4 
kbits/s,  not  shown  in  Fig.  7.6.2)  are  spare  ones. 

Data  acquisition  and  analysis  at  NORSAR 

A  data  acquisition  system  to  handle  the  data  from  the  two  stations  in 
Poland  hns  been  configured  at  NORSAR,  that  is  essentially  identical  to  the 
data  acquisition  systems  currently  handling  the  data  streams  from  the  regional 
arrays.  As  seen  in  Fig.  7.6.2,  the  communications  interface  is  NORSAR’s 
BUSC  system  (see  Paulsen  et  al,  1989).  Otherwise,  the  data  acquisition  system 
comprises  a  SUN  3/260  computer  with  a  Floating  Point  Accelerator,  a  1.2  Gb 
disk  drive  with  controller,  an  Exabyte  tape  drive  with  SCSI  controller,  and 
Ethernet.  This  system  will  make  it  possible  to  present  the  data  from  Poland 
to  the  IMS  in  the  same  way  data  from  the  arrays  are  made  available. 

The  24  October  1990  nuclear  test  at  Novaya  Zemlya  was  recorded  by  the 
new  seismic  station  at  Ksiaz,  and  the  recordings  for  both  the  short  period  and 
the  broad  band  stations  are  shown  in  Fig.  7.6.3. 


S.  Mykkeltveit 
R.  Paulsen 
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Fig.  7.6.1.  The  map  shows  the  location  of  the  regional  arrays  NORESS, 
ARCESS,  FINESA  and  GERESS,  as  well  as  the  two  new  stations  in  Poland. 
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Fig.  7.6.3.  Short  period  (sz,  sn,  se)  and  broad  band  (bz,  bn,  be)  3-component 
recording  at  Ksiaz  of  the  24  October  1990  nuclear  explosion  at  Novaya  Zemlya. 
The  upper  frame  shows  20  minutes  of  data,  whereas  the  lower  frame  covers  a 
one-minute  interval  around  the  P  onset. 
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7.7  Optimal  group  filtering  and  noise  attenuation  for  the 
NORESS  and  ARCESS  arrays 

A  generalization  of  Capon’s  (1970)  maximum  likelihood  technique  for  detec¬ 
tion  and  estimation  of  seismic  signals  has  recently  been  introduced  by  Kush- 
nir  et  al  (1990).  By  using  a  multidimensional  autoregressive  approximation  of 
seismic  array  noise,  a  technique  to  use  Capon’s  multichannel  filter  for  online 
processing  has  been  developed.  Such  autoregressive  adaptation  to  the  current 
noise  matrix  power  spectrum  has  been  shown  to  yield  good  suppression  of 
coherent  noise  processes.  This  paper  presents  further  results  of  our  studies  of 
of  adaptive  optimal  group  filtering  (AOGF)  which  now  has  been  applied  to 
more  than  80  earthquake  and  underground  nuclear  explosion  recordings  at  the 
small  aperture  arrays  NORESS  and  ARCESS  in  Norway  (Fig.  7.7.1). 

Compared  to  conventional  beamforming,  AOGF  is  demonstrated  to  pro¬ 
vide  high  signa!-to-noise  ratio  (SNR)  gains  (12-18  dB)  over  a  wide  frequency 
band  (0.2-5. 0  Hz),  while  retaining  an  undistorted  signal  waveform.  This  gain 
is  mainly  due  to  suppression  of  noise  at  low  frequencies,  where  both  noise 
power  and  noise  coherency  is  greatest. 

AOGF  has  also  been  tested  in  combination  with  statistically  optimal  algo¬ 
rithms  for  detection  and  onset  time  estimation  using  the  maximum  likelihood 
principle.  Here,  narrow-band  filtering  has  been  applied,  selecting  the  frequency 
band  most  appropriate  for  detecting  weak  arrivals  by  conventional  beamform¬ 
ing.  Even  in  this  case,  AOGF  provides  significant  gain  in  SNR  (typically  6-8 
dB  in  the  2-4  Hz  band)  compared  to  the  array  beam. 

Description  of  software  used 

The  conventional  technique  for  array  signal  processing  is  based  on  beam¬ 
forming,  band  pass  filtering  and  STA/LTA  detection.  Beamforming  is  the 
statistically  optimal  procedure  under  the  assumption  of  independent  noise 
among  the  M  different  receivers  of  the  array  and  in  this  case  it  improves 
power  signal-to-noise  ratio  by  a  factor  of  M.  For  correlated  noise  this  proce¬ 
dure  is  not  optimal  in  the  statistical  sense.  Bandpass  filtering  gives  additional 
improvement  of  SNR,  but  this  procedure,  evidently,  may  reduce  available  in¬ 
formation  about  the  signals.  The  STA/LTA  detector  is  a  simple  procedure, 
but  it  is  sensitive  only  to  variations  of  recorded  power  due  to  signal  arrival, 
and  it  does  not  make  use  of  spectral  differences  between  signal  and  noise.  So 
there  is  room  left  for  improvement  of  this  technique,  and  a  possibility  is  to  use 
statistically  optimal  procedures  accounting  for  noise  and  signal  features. 

Now,  at  NORSAR  a  technique  based  on  statistically  optimal  procedures 
has  been  tested.  This  technique  is  illustrated  in  Fig.  7.7.2  and  consists  of 
adapative  optimal  group  filtering  (AOGF)  followed  by  optimal  detection  (OD) 
and  optimal  signal  onset  time  estimation  (OE)  (see  Pisarenko  et  al,  1987). 
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AOGF  is  based  on  a  Wiener  optimal  filter  with  frequency  response 

H(  A)  =  G*(X)F~1(X)/G*{X)F~1(X)G(X) 

using  the  multichannel  array  record  Xt  =  St  +  &,  St  =  G*ut  for  input.  The 
output  of  this  filter  Yt  is  a  scalar.  AOGF  reduces  the  noise  power  to  a  mini¬ 
mum,  using  the  estimated  inverse  matrix  power  spectrum  F~1(X)  of  the  array 
noise  &.  The  vector  function  G(A)  represents  the  medium  frequency  response 
along  the  paths  from  the  seismic  source  to  the  array  sensors,  u*  is  the  scalar 
signal  of  the  seismic  source  (waveform).  For  a  homogeneous  medium,  is 
determined  by  signal  arrival  delays  at  different  stations.  If  AOGF  is  tuned 
properly  (G  -  G)  and  f  =  0,  then  the  AOGF  output  Yt  coincides  with  the 
signal  waveform  ut.  The  procedure  thus  avoids  introducing  any  distortion  of 
the  signal.  If  the  noise  is  uncorrelated  (F( X)  =  I)  (I  -  identity  matrix)  and 
G(A)  depends  on  delays  only,  then  AOGF  coincides  with  beamforming. 

The  inverse  power  spectrum  matrix  F~1(f)  is  estimated  from  pure  noise 
preceding  the  signal  by  a  procedure  of  adaptation  using  multichannel  autore¬ 
gressive  (AR)  modelling.  AR  estimation  of  this  large  size  matrix  function 
considerably  reduces  computations  and  has  experimentally  proved  to  be  an 
excellent  procedure  for  estimating  the  coherent  noise  matrix  power  spectrum. 

The  optimal  detector  (OD)  (Fig.  7.7.2)  is  based  on  autoregressive  mod¬ 
elling  and  moving  window  detection.  In  a  series  of  experiments  OD  proved  to 
exceed  STA/LTA  results  due  to  optimal  accounting  of  spectral  variations  in 
the  moving  window.  The  optimal  estimator  (OE)  (Fig.  7.7.2)  is  a  maximum 
likelihood  procedure  applied  to  the  time  interval  around  the  detection  time. 
It  is  likewise  based  on  an  autoregressive  modelling  procedure  and  is  sensitive 
to  small  variations  of  spectra  in  a  chosen  time  interval. 

Coherency  of  noise  and  optimal  group  filtering 

In  Kushnir  and  Lapshin  (1984)  it  is  shown  that  the  SNR  at  ihe  output  of 
OGF  tends  to  infinity  with  increasing  noise  coherency.  Coherent  noise  may 
be  represented  as  a  superposition  of  noise  contributions  from  several  spatially 
distributed  sources.  ARCESS  and  NORESS  noise  seems  to  have  a  strong 
coherent  component,  which  made  it  possible  to  achieve  power  SNR  gain  factors 
of  AOGF  relative  to  beamforming  of  about  70-80  at  ARCESS  and  20-30  at 
NORESS  (Kushnir  et  al ,  1990;  Kushnir  et  al ;  1989).  Coherent  noise  may  be 
recognized  from  spatial  spectral  FK-analysis,  where  noise  and  signals  coming 
from  different-  sources  correspond  to  different  peaks  in  the  slowness  vector 
plane.  Strong  coherent  noise  may  be  represented  in  the  frequency  domain  as 

N 

f(A)  =  £  £*(A )$k  +  erf(A),  A  e  (0,2;r] 

*=i 

where  &(A)  are  scalar  noise  source  signals,  <jf>(A)  is  the  vector  frequency  re¬ 
sponse  from  the  k- th  noise  source  to  the  array  receivers,  and  f?(A),  k  =  1, .  ,  N 
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are  independent  power  spectrum  matrices. 

m = E  aw&wSm  +  (i) 

,  k=l 

where  /*( A)  is  the  spectral  density  of  fo(A)  and  $(A)  is  the  matrix  power 
spectrum  of  diffuse  noise,  and  <j>*k  is  the  transposed  conjugated  vector.  It 
can  be  seen  that  if  e  -»•  0,  the  matrix  F( A)  projects  each  vector  signal  Q( A) 
recorded  by  the  array  to  the  “noise”  subspace,  i.e.,  linear  subspace 

{to,  k=  1 . N }: 

N 

/(A)Q(A)  £  CkMX) 

e-y  0  k=i 

and  ir_1(A)  projects  Q( A)  to  the  orthogonal  subspace  {<f>k(X)}k=:li  N.  Hence, 
the  power  of  each  signal,  coming  from  the  “noise”  subspace,  will  be  greatly 
reduced  by  OGF.  We  would  like  to  note  that  compensation  of  noise  by  OGF 
does  not  require  stationarity  of  the  noise  source  signals  £*(A),  but  stationurity 
of  <£(A)  is  important. 

Description  of  experiments 

Checking  of  noise  attenuation  stability.  This  paper  is  the  second  report 
concerning  our  attempts  to  use  AOGF  at  NORSAR.  In  the  first  one  (Kushnir 
et  al ,  1989),  we  described  theoretical  capabilities  of  AOGF  and  its  successful 
application  for  broad  band  extraction  of  one  of  the  smallest  nuclear  explosion 
signals  recorded  at  ARCESS.  But  in  the  case  described  in  Kushnir  et  al  (1989), 
adaptation  was  done  just  before  the  onset  time.  Therefore  some  doubts  re¬ 
mained  as  to  whether  AOGF  would  work  continuously  without  losing  its  good 
features  during  long  periods  of  filtering  time.  To  cheat  the  AOGF  performance 
stability,  we  applied  AOGF  to  some  long-term  records  of  pure  ARCESS  and 
NORESS  noise. 

In  Fig.  7.7.3  we  show  the  results  of  filtering  by  AOGF  the  25  channels 
of  a  40-  minute  long  NORESS  noise  record.  Only  the  first  2  minutes  of  the 
time  interval  were  used  for  adaptation.  AOGF  output  noise  power  is  gradually 
increasing,  but  SNR  power  gain  relative  to  the  beam  does  not  become  less  than 
20.  It  means  that  if  a  signal  arrives  during  this  time  interval,  its  SNR  will  be 
enhanced  by  AOGF  at  least  20  times  better  than  by  the  beam.  Looking  at 
Fig.  7.7.3,  one  may  predict  that  OGF  can  work  much  better  than  the  beam 
during  one  hour  or  even  more. 

We  applied  AOGF  to  5  ARCESS  and  5  NORESS  records  each  consisting 
of  25  minutes  of  pure  noise,  computing  the  SNR  power  for  each  5  minutes. 
The  results  are  presented  in  Fig.  7.7.4.  From  Fig.  7.7.4b  we  see  that  the  gain 
of  AOGF  relative  to  beamforming  for  NORESS  is  between  15  and  25,  for 
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ARCESS  from  30  to  75.  Again,  adaptation  is  done  during  the  first  2  minutes, 
and  that  is  why  we  have  the  largest  gain  for  the  first  5  minutes,  including  the 
interval  of  adaptation.  AOGF  includes  two  stages:  adaptation  and  filtering. 
The  first  is  time  consuming,  the  second  is  not.  If  readaptation  is  only  needed 
once  an  hour  or  so,  one  may  attempt  to  use  several  AOGFs,  tuned  to  different 
azimuths  and  apparent  velocities  in  the  same  manner  as  done  by  beamforming 
at  NORSAR,  and  this  may  be  more  effective  for  noise  suppression  in  online 
processing.  An  illustration  of  such  an  application  of  AOGF  is  shown  in  Fig. 
7.7.5,  where  the  gain  of  AOGF  is  presented  versus  azimuth.  The  velocity  is 
constant  and  equal  to  13  km/sec.  The  peak  of  the  gain  corresponds  to  the 
southern  direction. 

Extraction  of  missed  phases.  The  capability  of  AOGF  to  work  for  a  long  ti.  .e 
without  readaptation  also  makes  it  possible,  using  pure  noise  before  P  arrivals 
for  adaptation,  to  enhance  the  SNR  of  seismic  phases  located  in  time  far  from 
the  P  arrival.  Two  examples,  described  below,  illustrate  the  application  of  all 
three  optimal  algorithms:  adaptive  optimal  group  filtering,  optimal  detection, 
and  onset  time  estimation  for  extraction  of  missed  seismic  phases.  It  some¬ 
times  happens  that  regional  phases  like  Pn  or  Lg  are  missed  by  the  detector 
applied  at  NORSAR  in  the  regular  processing  of  NORESS  and  ARCESS  data. 
We  tried  to  use  AOGF,  OD  and  OE  to  remedy  this  situation  for  some  of  these 
cases,  and  they  seemed  to  work  well.  Two  cases  were  considered. 

a)  For  case  1  (Fig.  7.7.6a),  the  Lg  phase  was  detected  on  NORESS,  but 
the  Pn  and  Sn  phases  were  missed.  Fig.  7.7.6b  shows  data  for  the  same 
event,  as  recorded  on  the  FINESA  array  in  Finland.  In  Fig.  7.7.6c  the 
result  of  AOGF,  OD  and  OE  and  beamforming  are  shown  for  a  window 
containing  the  presumed  Pn  onset.  The  onset  time  of  the  Pn  wave  is 
clearly  seen,  and  the  maximum  of  the  OE  curve  is  close  to  the  theoretical 
Pn  arrival.  (The  theoretical  Pn  arrival  time  is  computed  from  the  event 
location  derived  from  the  FINESA  array  data.)  The  theoretical  and 
experimental  P  arrivals  are  shown  on  Fig.  7.7.6a.  The  first  detection 
shown  in  Fig.  7.7.6d  is  questionable  (is  it  caused  by  the  Sn  phase  from  the 
event  considered  or  not)  because  the  theoretical  Sn  arrival  time  differs 
from  the  experimental  one.  The  second  detection  in  Fig.  7.7.6d  is  more 
likely  related  to  the  Sn  phase  from  the  event  considered. 

b)  For  case  2  (Fig.  7,7.7a)  the  regular  NORESS  processing  detected  the  Pn 
and  Sn  phases,  but  no  Lg  was  detected.  This  interpretation  is  confirmed 
by  inspecting  the  ARCESS  records.  We  have  done  an  adaptation  for 
AOGF  just  before  the  Pn  arrival.  The  results  are  shown  in  Fig.  7.7.7b, 
where  all  three  phases  Pn,  Sn  and  Lg  are  evidently  detected.  Their 
arrival  times  are  shown  on  Fig.  7.7.7a  as  well. 

These  results  suggest  possible  applications  of  AOGF  in  combination  with 
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OD  and  OE  for  improvement  of  the  detection  lists  generated  by  NORSAR’s 
regular  processing  of  array  data. 

Adaptive  group  filtering  of  seismic  signals 

AOGF  was  applied  to  the  processing  of  35  teleseismic  earthquakes  and 
44  nuclear  explosion  signals  recorded  at  NORESS.  Earthquake  locations  are 
shown  in  Fig.  7.7.8.  Most  of  the  explosions  are  located  at  the  Shagan  River 
site.  This  selection  of  events  was  made  with  a  view  to  future  work  on  the 
discrimination  problem.  In  all  the  cases,  the  SNR  was  much  improved  by 
AOGF,  but  the  events  may  be  divided  into  two  classes.  One  class  is  where 
SNR  was  small  in  the  first  place  and  the  other  where  it  was  large.  In  the  first 
class  (example  in  Fig.  7.7.9)  improvement  of  SNR  by  AOGF  is  evident,  and 
in  the  other  example  (Fig.  7.7.10a)  the  improvement  is  not  evident  in  a  broad 
frequency  band  (0.2-5.0  Hz).  But  if  we  consider  only  the  low  frequency  band 
(for  example,  0.2-0.5  Hz)  the  advantage  of  the  group  filter  is  seen  clearly  (Fig. 
7.7.10b).  This  is  due  to  two  factors:  The  first  is  that  at  this  latter  frequency 
band  the  SNR  is  worse  for  the  short-period  records  used,  and  the  second  is 
that  the  gain  of  AOGF  seems  to  be  largest  due  to  the  high  noise  coherency  in 
this  band. 

Conclusion 

By  application  of  advanced  techniques  (apative  optimal  group  filtering,  op¬ 
timal  detector,  optimal  estimator),  based  on  statistical  algorithms,  processing 
of  NORESS  and  ARCESS  array  data  may  be  improved.  This  improvement 
could  be  achieved  by  taking  optimal  advantage  of  array  noise  coherency  and 
signal  frequency  contents.  Compared  to  conventional  beamforming,  AOGF  is 
demonstrated  at  NORESS  and  ARCESS  to  provide  SNR  gains  from  12  to  18 
dB  over  a  wide  frequency  band  from  0.2  to  5  Hz,  while  retaining  an  undis¬ 
torted  signal  waveform.  This  gain  is  mainly  due  to  suppression  of  noise  at  low 
frequencies,  where  both  noise  power  and  noise  coherency  is  greatest.  But  even 
for  the  frequency  band  from  2  to  4  Hz,  AOGF  provides  a  typical  gain  in  SNR 
of  the  order  6  to  8  dB  compared  to  the  array  beam.  The  performance  of  AOGF 
appeared  to  be  stable  in  time  without  need  for  readaptation,  for  time  inter¬ 
vals  at  least  one  hour  long.  This  may  result  in  possible  application  of  AOGF 
to  supplement  conventional  beamforming  for  online  signal  detection  and  ex¬ 
traction.  The  proposed  technique  may  be  utilized  also  in  post-processing  for 
detecting  and  onset  time  estimating  of  seismic  phases  missed  by  conventional 
methods.  As  a  result,  it  may  reduce  the  detection  threshold  and  improve  the 
reliability  of  event  detection  and  location  using  array  signals. 
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Fig.  7.7.1.  The  locations  and  array  geometry  for  the  NORE3S  and  ARCESS 
arrays.  The  two  arrays  are  essentially  identical,  each  having  24  elements  in 
four  concentric  rings  (A,  B,  C,  D)  plus  a  center  element  (AO).  The  diameter 
of  the  outer  (D)  ring  is  3.0  km.  There  are  3-component  seismometers  at  AO 
and  three  of  the  seven  sites  in  the  C-ring. 
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Fig.  7.7.2.  Statistically  optimal  algorithms  for  preliminary  array  signal  pro¬ 
cessing.  X*  —  multichannel  seismic  record;  St  —  signal;  —  noise;  H  =  H(X) 
—  frequency  response  of  OGF;  F_1  =  F_1(A)  —  inverse  matrix  power  spec¬ 
trum;  G  =  G(A)  —  vector  medium  frequency  response;  yt  —  scalar  output 
of  AOGF,  i?(i)  —  detector  statistic;  y/+m  —  observations  in  moving  win¬ 
dow  (width  =  m  samples);  A  —  asymptotically  sufficient  statistics;  C~l  — 
inverse  covariance  matrix  of  A  computed  at  adaptation  step;  Q{1)  —  like¬ 
lihood  function  of  onset  time  fe;  LF  —  likelihood  functions  of  observations 

Y{  =  (Y, . Y,),r#t  =  (y<+, . Ym). 
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Gain  of  optimal  group  filter 


Fig.  7.7.4*  Power  SNR  gain  of  optimal  group  filter  computed  each  five  min¬ 
utes:  a)  relative  to  beamforming;  and  b)  relative  to  average  NORESS  channel 
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Fig.  7.7.5.  Group  filter  gain  for  pure  noise  versus  azimuth.  The  apparent 
velocity  is  equal  to  13  km/sec.  The  maximum  corresponds  to  the  southern 
direction. 
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NORESS  data.  No  Pn  detection  on  NORESS.  Only  Lg  detected. 
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Fig.  7.7.0a.  Illustration  of  the  AOGF,  OD  and  OE  procedures  as  applied  to 
NORESS  data  for  a  small  event  in  Estonia,  from  which  only  the  Lg  phase  was 
detected  by  the  regular  processing  at  NORSAR  of  NORESS  data.  The  arrows 
show  theoretical  and  computed  (by  the  OE  procedure)  arrival  times  for  Pn 
and  Sn  phases.  Circles  mark  the  time  interval,  in  which  there  was  a  search  for 
a  Pn  phase  (see  Fig.  7.7.6c),  triangles  mark  the  interval  for  searching  for  the 
Sn  phase  (see  Fig.  7.7.6d). 
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Fig.  7.7.8c.  NORESS  Pn-wave  detection  and  onset  estimation  performed 
on  AOGF  output  data.  “V”  and  “AZI”  coi respond  to  apparent  velocity  and 
azimuth  used  in  AOGF,  “Parr”  is  computed  by  OE  relative  to  the  beginning 
of  the  interval  for  the  search  for  the  Pn  phase. 


NORESS  data.  No  Lg  detected  on  NORESS. 
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Fig.  7.7.7a.  Illustration  of  the  AOGF,  OD  and  OE  procedures  as  applied  to 
an  event  from  the  USSR/Finland  border  region,  for  which  the  Lg  phase  was  not 
detected  by  the  NORESS  regular  processing.  (The  Sn  phase  at  approximately 
09.04  was  detected,  but  incorrectly  assigned  as  Lg.)  Arrows  show  estimates  of 
arrival  times,  using  the  OE  procedure. 
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Fig.  7.7.8.  Map  showing  the  location  of  35  earthquakes  used  in  this  study. 
Most  explosions  used  are  located  at  the  Shagan  River  site. 
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Fig.  7.7.10.  Example  of  nuclear  explosion  record  filtered  by  beam  and  AOGF 
in:  a)  broad  frequency  band  and  b)  low  frequency  band. 


